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LES PROPOS DU SECRETAIRE 


1) Notre bulletin n’était jusqu’a présent qu’un tiré a part du Bulletin 
d’information de ’ UGGI. 

Nous avions accepté cette fagon de faire, espérant pouvoir réduire les 
frais de publication. L’expérience a fait apparaitre certains inconvénients 
du systéme et notamment la difficulté de publier dans un bulletin s’adressant 
a l'Union de multiples et longs exposés dont |’intérét était limité aux membres 
de l’Association d’Hydrologie Scientifique. 

Le secrétaire général de l’Union, Mr Laclavére, a trés bien compris ces 
difficultés et c’est avec son accord que nous aurons dés 1957 un bulletin indé- 
pendant dont ces lignes préfacent le premier numéro. 

Comme on I|’a décidé 4 Rome, ce bulletin présentera une partie adminis- 
trative et une autre réservée aux publications scientifiques. Les charges finan- 
ciéres nous obligeront sans doute parfois a limiter la partie scientifique : nous 
ferons de notre mieux a cet égard. Nous n’avons pas 4a revenir sur les raisons 
de la nécessité de la partie dite administrative. Elle serait justifiée rien que 
par la nécessité de vous tenir au courant de la vie des Associations et Comités 
qui ont des rapports avec nous, tout en n’appartenant pas a notre Union. 
L’organisation de nos Assemblées et des Symposia sera d’ailleurs facilitée 
par les communications que nous vous adresserons. 


2) Dans ces propos du secrétaire, vous nous permettrez lors de la parution 
de chaque numéro, de vous présenter en quelques mots, des communications, 
des avis, des notes de tous genres qu’il est impossible de mettre sous un titre 
autre que celui de « propos ». Nous vous demanderons d’y préter quelque 
attention si vous voulez vous tenir au courant de la vie de 1’Association. 
Nous nous permettrons au cours de ces propos, de vous parler assez familiere- 


ment, comme dans une conversation. 


3) Nous profitons'immédiatement de cette liberté que nous nous accor- 
dons pour vous demander de nous conseiller sur le point suivant: 1’Union, 
dans sa politique actuelle, pousse a l’organisation par les Associations de 
Symposia, dans le genre des Symposia Darcy 4 Dijon, en dehors des Assemblées 
Générales. Ces Symposia attirent l’attention sur nos travaux dans des contrées 
qui ne connaissent guére nos Assemblées par suite de la limitation du nombre 
de celles-ci. 

L’UNESCO se montre d’ailleurs favorable a l’organisation de ces Sym- 
posia, mais le secrétaire-trésorier de l1’Association voit ses finances menacées 
par le coat de la publication des mémoires présentés a ces Symposia. La 
générosité de 1’ UNESCO s’attache surtout a favoriser les déplacements des 
participants. Or, l’organisation et les publications du Symposia de Dijon 
nous coiiteront environ 6000 dollars ! Le Conseil, réuni 4 Dijon, a bien voulu 
accepter notre proposition de ne faire aucune distribution gratuite de ces 
« Transactions », mais en dépit des efforts de certains adhérents dont je vou- 
drais parler ci-dessous, l’opération « Dijon » restera vraisemblablement deéfi- 
citaire au point de vue financier, bien que les lettres que nous avons recues de 
partout en font une opération des plus heureuses au point de vue scientifique. 
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4) Je viens de dire un mot de ceux d’entre vous qui se sont dévoués pour 
assurer la diffusion des Comptes-Rendus des Symposia de Dijon. Nous ne 
pouvons nous permettre de citer des noms, mais nous adressons les remer- 
ciements de l’Association a tous ceux qui ont fait connaitre ces publications, 
notamment par des articles appropriés dans leurs revues nationales. 

A Rome et a Dijon, nous avons maintes fois répéte: on ne connait pas 
suffisamment notre travail. Nous vous demandons encore une fois votre aide 
4 tous. Nous avons l’impression que de courts comptes-rendus dans les jour- 
naux scientifiques nationaux ont un rendement des plus favorables sur la vente. 

Nous nous excusons de ces considérations bien terre a terre: nous ne 
pensions pas avoir une 4me de commergant, mais nous apercevons plus que 
jamais que nous ne pouvons vivre qu’avec les deux pieds sur terre. 


ASSEMBLEE DE TORONTO 


1) Nous avions fixé le 1°" avril comme date limite de la rentrée des titres 
des communications qui seront présentées. Nous l’avons fait par avis per- 
sonnel aux représentants nationaux qui nous ont été désignés. Nous l’avions 
rappelé dans le numéro du 1¢* Décembre 1956 du Bulletin. Mais au moment 
ou nous rédigeons cet avis le numéro en question n’est pas encore sorti de 
presse. 

Nous craignons donc que certains pays aient oublié la date limite. Nous 
accepterons donc les relevés de titres jusqu’au 30 avril prochain. 


2) Nous demandons au représentant pour lhydrologie de chaque pays 
de ne pas oublier la date du 1°" juin 1957 pour la rentrée des textes des rapports. 
Pouvons-nous insister pour que les auteurs n’oublient pas de joindre a leur 
rapport un résumé comportant au maximum une page dactylographiée avec 
double interligne. Si c’est possible, le représentant national nous enverra ces 
résumés aussit6t que possible, méme avant le texte des rapports. 


3) A titre purement indicatif, nous donnons ci-dessous un premier essai 
du programme, en insistant beaucoup sur le mot essai. 


Mardi 3 septembre 


Avant-midi : Assemblée Générale de |’?Union. Convocation Hall. 
Apres-midi: Séance du Conseil de |’Association. 


Mercredi 4 septembre 


Avant-midi: Séance commune avec les Associations d’Océanographie et de 
Metéorologie : Théme général: Bilan Hydrologique. 
L’Association d’Hydrologie interviendra en présentant celles des 
communications se rattachant a ce theme et vraisemblablement 


__ celles sur influence de la végétation. 
Apres-midi: Idem. 


Jeudi 5 septembre 


Avant-midi: Com. Neiges et Glaces: 1" Séance. 
ae Com. Eaux souterraines: 1° Séance. 
Apres-midi : Assemblée générale. Adresse présidentielle. Rapport du Secré- 
taire. Exposé de diverses questions et discussions. 
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Vendredi 6 septembre 


Avant-midi: Com. Neiges et Glaces: 2™e Séance. 
Com. Eaux souterraines : 2™e Séance. 
Aprés-midi: Com. Neiges et Glaces: 3™¢ Séance. 
Com. Eaux souterraines: 3™° Séance. 


Samedi 7 septembre 


Avant-midi: Com. Neiges et Glaces: 4™e Séance. 
Com. Eaux souterraines: 4™¢ Séance. 


Dimanche 8 septembre 


Excursion. 


Lundi 9 septembre 


Avant-midi: Com. Neiges et glaces: 5™® Séance. 
Com. Eaux souterraines: 5™e Séance. 

Apres-midi: Com. Erosion Continentale : 1"¢ Séance. 
Comité des Mesures et Instruments. 


Mardi 10 septembre 


Avant-midi: Com. Erosion Continentale : 2™® Séance 
Comité des Précipitations: 17° Séance 
Apres-midi: Com. Erosion Continentale : 3™¢ Séance 
Com. des Eaux de Surface: 1° Séance (Evaporation) 


Mercredi 11 septembre 


Avant-midi: Com. Eaux de surface: 2™¢ Séance (Evaporation) 
Comité des Neiges et Glaces: 6™® Séance 
Aprés-midi: Com. Eaux de surface : 3™° Séance 
Comité des Neiges et Glaces: 7™° Séance 


Jeudi 12 septembre 


Avant-midi: Rosée, condensations et précipitations occultes (1"¢ Séance) 
Aprés-midi: Rosée, condensations et précipitations occultes (2™¢ Séance) 
Couwité des précipitations : 2° Séance 


Vendredi 13 septembre 


Avant-midi: 9 h. Com. des Eaux de surface (4™° Séance) 
a 11 h. Discussion sur la Standardisation des caractéristique 
hydrologiques. 

Aprés-midi; 14.30 h. Com. Eaux de surface : 5™° Séance. 
17 h. Assemblée générale. 


Samedi 14 septembre 


Avant-midi: Assemblée de cloture de I’Union. 


. THE SECRETARY’S PROPOSALS 


1) Up to the present our Bulletin has been printed only as part of the 
Information Bulletin of UGGI. We have accepted this method in the hope 
that publication costs would be reduced. Experience has shewn certain 
difficulties particularly the publication in the Union’s Bulletin of a large number 
of papers of limited interest to members of the Association of Scientific 
Hydrology. 

The Union’s General Secretary, Colonal Laclavére, has been aware of 
our difficulties and has agreed to our proposal to publish independently in 
1957 along the following lines. 

In accordance with a decision made in Rome, this Bulletin will consist 
of two parts, one administrative and the other reserved for scientific 
publications. For financial reasons, we may be compelled at times to limit 
the scientific part, but we shall avoid this if possible. It is unnecessary to 
reiterate the reasons for the administrative part. This would be justified only 
in keeping you informed of the activities of other Associations and Committees, 
which are not members of our Union. The organisation of our Meetings and 
Symposia will be made easier by these communications. 


2) The Secretary proposes to present in each number short com- 
munications, notices, etc., under the general heading of « Remarks ». I would 
ask you to pay particular attention to these if you wish to be kept informed 
of the Association’s affairs. In the course of these remarks I shall take the 
liberty of adressing you informally. 


3) I would ask you to make note of the following points. The present 
policy of the Union is to extend the organisation of Symposia by the 
Association, as for example the « Darcy » Symposia at Dijon, outside the 
General Assemblies. These Symposia draw attention to our work in those 
countries where our Assemblies are scarcely known, due to their infrequency. 

UNESCO favours the organisation of these Symposia, but the Asso- 
ciation’s Treasurer is concerned about the cost of publication of these 
proceedings. The assistance of UNESCO is directed to the publication of 
these reports. The organisation and publications of the Dijon Symposia will 
cost about 6,000 dollars. At Dijon, the Council agreed to the proposal that 
the Transactions should not be distributed without payment. In spite of the 
efforts of some members, about whom I would like to say more later, the 
Dijon meeting will probably result in a financial loss. On the other hand, 


letters which I have received indicate that it was a most successful meeting 
from a scientific point of view. 


4) I have just referred to those who have assisted in the distribution of 
the reports of the Dijon Symposia. I cannot mention individuals by name, 
but the thanks of the Association are due to those members who have made 
known these publications, particularly in their national journals. 

In both Rome and Dijon, it was repeatedly said that our work is 
insufficiently well known. I ask you all again for your help. I believe that 
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short articles in national scientific journals have a favourable effect on our 
sales. 

I apologise for these mundane. considerations. I do not wish to raise 
commercial considerations, but our Association can only thrive if our feet 
are planted firmly on the ground. sr 


TORONTO ASSEMBLY 


1) April Ist was fixed as the last date for the return of titles of papers. 
National representatives were advised of this personally. I drew attention to 
this in the December 1956 Bulletin, but at the time of drafting this notice the 
total number of papers is not available. It would appear that certain countries 
have overlooked this date, and I shall accept titles until 30th April. 


2) National representatives are asked to let me have their reports not 
later than Ist June, 1957, an I would ask that authors should not forget to 
include with their Reports a summary of a single page, typewritten and double 
spaced. These summaries should be forwarded by national representatives 
before the reports, if possible. 


3) As a guide, a first draft of the programme of the various Commissions 
is given, subject to alteration. 


Tuesday 3rd September 
Morning: General Meeting of the Union. Convocation Hall. 
Afternoon: Council Meeting of the Association. 


Wednesday 4th September 
Morning: Ordinary Meeting with the Oceanographic and Meteorological 
: Associations. General theme: Hydrological balance. 
The Hydrological Association will present some paper on this 
subject and probably on the influence of vegetation. 
Afternoon: Ditto. 


Thursday 5th September 

Morning: Snow and Ice: Ist session. 
Ground Water : Ist session. 

Afternoon: General Meeting. Presidential address: Secretary’s Report 
Outline of various questions and discussions. 


Friday 6th September 

Morning: Snow and Ice: 2nd session. 
Ground Water: 2nd session. 

Afternoon: Snow and Ice: 3rd session. 
Ground Water: 3rd session. 


Saturday 7th September 
Morning: Snow and Ice: 4th session. 
Ground Water: 4th session. 


Sunday 8th September 
Excursion. 


Monday 9th September 
Morning: Snow and Ice: Sth session. 
Ground Water: 5th session. 
Afternoon: Continental Erosion: 1st session. 
% Measures and Instruments. 


Tuesday 10th September 

Morning: Continental Erosion: 2nd session. 
Rainfall (1st session). 

Afternoon: Continental Erosion: 3rd session. 
Surface Water: Ist session (Evaporation). 


Wednesday 11th September 

Morning: Surface Water: 2nd session (Evaporation). 
Snow and Ice: 6th session. 

Afternoon: Surface Water: 3rd session. 
Snow and Ice: 7th session. 


Thursday 12th September | 


Morning: Dew, condensation and similar precipitations: Ist session. 


Afternoon: Ditto: 
Rainfal (2nd session). 


Friday 13th September 
Morning: 9 a.m. Surface Water: 4th session. 


11 a.m. Discussion on standardisation of hydrological charac- 


teristics. 
Afternoon: 2.30 p.m. Surface Water: 5th session. 
5 p.m. General Meeting. 


Saterday 14th September 
Morning: Closing Meeting of the Union. 


UNESCO/NS/140 
PARIS, le 17 janvier 1957 
Traduit de l’anglais 


ORGANISATION DES NATIONS UNIES 
POUR L’EDUCATION, LA SCIENCE ET LA CULTURE 


PROGRAMME DE LA ZONE ARIDE 


RAPPORT SUR LA ONZIEME SESSION 
DU COMITE CONSULTATIF DE RECHERCHES SUR LA ZONE ARIDE 


University College (Canberra) 
13-16 octobre 1956 
et 
Melbourne 
25 octobre 1956 


1. Le Comité consultatif de recherches sur la zone aride a tenu sa onzieéme 
session du 13 au 16 octobre 1956 a l’University College de Canberra, et le 
28 octobre a Melbourne. 


2. L’Organisation des Nations Unies, les Institutions spécialisées et les 
organisations internationales non gouvernementales qui s’intéressent aux 
recherches scientifiques et techniques sur la zone aride avaient été invitées 
a envoyer des représentants. Le Directeur général a également invité les membres 
australiens des colléges de consultants spécialisés dans les recherches sur la 
zone aride a participer a la session, qui était publique. La liste des participants 
est donnée en annexe au présent rapport. 


3. La session a été ouverte par M.G. Aubert (France), Président de la 
dixiéme session, qui a remercié le Gouvernement australien d’avoir invité 
lV’Unesco a tenir cette session a Canberra, et le Président et le Conseil de 
University College d’avoir mis leurs locaux a la disposition du Comite. 


4. Le Comité a porté M. H. O’R. Sternberg (Brésil) a la présidence pour 
la durée de la onziéme session. Le Président a constitué un sous-comité (com- 
posé de MM. G. Aubert (France), M.H. Ganji (Iran) et Herbert Greene 
(Royaume-Uni) qu’il a chargé de rédiger le rapport sur la session. 


5. Le Comité a souhaité la bienvenue a M. C. W. Thornthwaite (Etats- 
Unis d’Amérique) qui lui a été adjoint depuis la derniére session, ainsi qu’a 
M. M.H. Ganji, suppléant de M. A. Behnia (Iran), empéché d’assister a 
la session. Il a regretté que le troisiéme de ses nouveaux membres, M. G.V. Bo- 
gomoloy (Union des Républiqués soviétiques socialistes), venant de Moscou, 
ait été retardé en cours de route. Le Comité a d’autant plus déploré l’absence 
d’un représentant de |’Organisation des Nations Unies que le Conseil écono- 
mique et social a, au cours de sa vingt et uniéme session (avril-mai 1956), 
adopté des résolutions chargeant le Secrétaire général d’entreprendre une action 
dans certains domaines auxquels le Comité porte un vif intérét. 


6. Le Comité a pris note avec une grande satisfaction des dispositions 
que le Directeur général et le Gouvernement australien ont adoptées en vue 


_. d@’organiser 4 Canberra, du 17 au 20 octobre 1956, un colloque sur la clima- 


tologie de la zone aride, notamment sur la micro-climatologie, et des mesures 
que le Gouvernement australien a décidées pour permettre aux participants 
a ce colloque de faire ensuite un voyage d’étude. Le Comité a prié le Directeur 


9 


t 


général d’exprimer sa gratitude au Gouvernement australien, a la Common- 
wealth Scientific and Industrial Research Organization (que le Gouvernement 
avait chargé de préparer le colloque), a M. B. T. Dickson, Président, et aux 
membres du comité local d’organisation. 


7. Le Comité a approuvé le rapport du Secrétariat, qui rend compte de 
oeuvre accomplie dans le cadre du programme de la zone aride depuis Ja 
session précédente. II s’est déclaré satisfait de l’action exercée par le Secrétariat 
au cours de cette période. Il a remercié le représentant de l’Organisation 
météorologique mondiale qui lui a transmis les voeux de son Président, en 
réaffirmant l’intérét que l’}OMM porte aux travaux du Comité et le désir 
qu’elle a d’y prendre part. Il a également remercié de ses voeux le Conseil 
scientifique pour |’Afrique au sud du Sahara et a pris note avec satisfaction 
de l’activité exercée par ce Conseil dans le domaine des recherches sur la zone 
aride. Il a exprimé l’espoir que le Directeur général invitera désormais ce Con- 
seil a se faire représenter aux sessions du Comite. 


Projet majeur 


8. Le Comité consultatif a étudié avec soin et dans le détail les propositions 
concernant un projet majeur relatif aux recherches scientifiques sur les terres 
arides, que le Directeur général a présentées dans le Projet de programme 
et de budget de l’Unesco pour 1957-1958, dont sera saisie la Conférence géné- 
rale de 1’Unesco lors de sa neuviéme session, 4 New Delhi, en novembre 1956. 
Le Comité a formulé les recommandations suivantes : 


Généralités 
Recommandation | 


9. Il ressort nettement des débats de la Conférence des Nations Unies 
sur la conservation et l’utilisaticn des ressources naturelles (1948) que les 
problémes des régions arides et semi-arides ont été négligés. C’est pourquoi 


l’Unesco a accepté la mission d’orienter |’attention des milieux scientifiques 
vers ces questions. 


10. Les zones arides et semi-arides s’étendent sur de vastes territoires. 
Les conditions physiques y grévent lourdement la vie végétale et animale. 
Les groupes humains, qui luttent pour survivre en milicu hostile, y sont faibles 
et isolés. Ils n’ont probablement pas assez — ou pas du tout —— d’>hommes 
de sciences et pourtant il n’est pas possible de résoudre leurs problémes — 
ni méme de les rendre moins aigus — autrement que par un effort continu 
faisant appel a la collaboration de plusieurs sciences. A ce prix, on devrait 
pouvoir accomplir des progrés sur un vaste front. 


11. Le comité consultatif de recherches sur la zone aride, considérant la 
modicité des ressources dont dispose l’Unesco, a rejeté les propositions tendant 
a les consacrer a la création d’un institut ou d’un laboratoire international 
de recherches sur la zone aride, qui risquerait d’en absorber |’intégralité. 
Le Comité consultatif a décidé que ce qu’il peut faire de plus utile est de: 

a) Rassembler et diffuser des informations en dressant l’inventaire des 
recherches et en organisant des colloques; 

b) Faciliter la planification, la coordination et l’exécution des recherches 
sur la zone aride par l’intermédiaire d’instituticns reconnues et grace au 
concours de personnes spécialement choisies: 


c) Accorder des bourses pour la formation de spécialistes en matiére de 
recherches sur la zone aride. 
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12. Une importante addition au programme a été approuvée a l’issue. 
d°’une réunion tenue a Socorro (avril-mai 1955) conjointement avec |’ American 
Association for the Advancement of Science. Il a été observé que, méme 
dans les pays oti la valeur des ressources scientifiques est amplement reconnue, 
on néglige incontestablement d’appliquer aux problémes des terres arides 
les connaissances acquises. Ce fait est imputé non seulement au manque de 
coordination entre les diverses sciences, mais aussi au manque de contact 
entre les savants et les profanes. Il a été décidé qu’il y avait lieu de s’employer 


* davantage a créer des comités locaux et nationaux comprenant de nombreuses 


personnalités non scientifiques. 


13. Le Comité consultatif est convaincu qu’un grand pas sera fait dans 
cette direction si les Postes régionaux de coopération scientifique qui doivent 
faciliter la réalisation du projet majeur, considérent qu’ils ont pour mission 
primordiale d’encourager et d’appuyer |’action de ces comités locaux et natio- 


X 


naux. Des crédits additionnels seront nécessaires a cet effet. 


14. Les membres des Postes de coopération scientifique, en nouant des 
relations cordiales avec les membres des comités locaux et nationaux, doivent 
provoquer des demandes raisonnables d’assistance en vue de projets de recher- 
ches s’inscrivant dans le cadre du programme relatif a la zone aride. Des pro- 
jets bien concus, émanant des comités locaux ou nationaux, et aboutissant 
a des résultats tangibles, peuvent trés efficacement contribuer 4 éliminer les 
obstacles aux progrés. Une fois bien comprise des milieux locaux, la tache 
entreprise par V’'Unescd pourrait étre puissamment soutenue par l’une ou 
Tautre des trés riches collectivités installées en zone aride. 


15. Les conditions ainsi créées faciliteraient l’établissement d’un prce- 
gramme judicieusement coordonné dans lequel les comités locaux et natio- 
naux d’une part, et 17?Unesco, d’autre part, se répartiraient la tache. 


Recommandation 2 


16. Le Comité, considérant que les limites des régions arides et semi- 
arides sont dynamiques plut6t que statiques, que de nouveaux Etats pourront 
devenir membres de |’Unesco au cours des cinq ou six prochaines années, 
et que, dans la liste incorporée au préambule du projet majeur plusieurs 
territoires non autonomes ne sont pas mentionnés, a recommandeé d’éliminer 
du texte du projet majeur la liste nominative des pays. 


A. — Sessions du Comité consultatif 


Recommandation 3 


17. Considérant qu’il importe d’entreprendre aussit6t que possible en 
1957, la tache prévue en faveur du Moyen-Orient dans le cadre du projet 
majeur, et tenant compte de la proposition tendant a ce que |’exécution de 
ce projet se poursuive pendant cing ou six ans, le Comité a recommandé que 
le Directeur général envisage a titre exceptionnel de convoquer deux sessions 
du Comité en 1957. 


B — Comités nationaux et locaux de coopération 


Recommandation 4 


18. Le Comité a recommandé que des mesures soient prises pour faire 
en sorte que des exemplaires du Guide des travaux de recherche sur la mise 
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en valeur des régions arides (Guidebook to research data for arid zone develop- 
ment) soient mis a la disposition des comités nationaux et locaux de coopé- 
ration, des instituts de recherches sur la zone aride, ainsi que d’organisations 
et de personnalités spécialement choisies qui pourraient assurer |’utilisation 
rationnelle de cet ouvrage. 


Recommandation 5 


19. Le Comité a recommandé que le Secrétariat publie un bulletin pério- 
dique ayant pour principal objet d’établir un lien entre les Comités de coopeé- 
ration et le Comité consultatif. Ce bulletin contiendrait des renseignements 
sur les activités entreprises dans le cadre du programme, sur les réunions 
scientifiques dont les travaux touchent a la zone aride, notamment celles 
auxquelles participent des membres du Comité consultatif, et donnerait des 
informations sur les comités de coopération. I] fournirait aussi des indications 
sur les bourses. Le premier numéro contiendrait un article exposant la fagon 
dont se constitue un comité national ou local de coopération. Le montant 
nécessaire a la publication de ce bulletin pourrait étre prélevé sur les crédits 
inscrits au budget aux fins de publication d’une deuxi¢me édition du 
Directory of Institutions engaged in Arid Zone Research (Répertoire des insti- 
tutions s’occupant de recherches sur la zone aride) car, de l’avis du Comité, 
ce répertoire serait peu demandé, a en juger par le faible chiffre de vente de 
la premiére édition. 


C — Rapports-inventaires sur les recherches 


Recommandation 6 


20. Le Comité a recommandé qu’en 1957 un contrat soit passé pour 
la rédaction d’un ouvrage retracant l’historique du mode d’utilisation des 
terres et notamment de l’agriculture, dans les régions arides et semi-arides : 
cet ouvrage prendrait place dans la collection des rapports-inventaires non 
périodiques. Le Comité attache une grande importance a ce travail. 


Recommandation 7 


_ 21. Le Comite a recommande qu’un contrat soit passé en 1957 pour rédac- 
tion, aux fins de publication dans la collection des inventaires périodiques, 
de deux publications dont la premiére serait un inventaire critique sélectif 
des principales contributions apportées, depuis le Colloque d’Ankara de 1952, 
au domaine de l’hydrologie, du point de vue particulier de la zone aride, et 
la seconde un inventaire analogue des travaux effectués au sujet des plantes 
médicinales de la zone aride. Le Comité a pris note avec satisfaction de l’offre 
qu’a faite l’Organisation mondiale de la santé d’aider le Comité a établir 
ce dernier inventaire. Ces deux inventaires, qui ne comprendraient qu’une 
quarantaine de pages imprimées, seraient publiés séparément. 


D — Colloques 


Recommandation 8 


PT, Alte Comité a pris note avec un vif intérét du projet d’organisation, 
au cours du troisieme trimestre 1957, d’un colloque sur la zone aride sous 
les auspices du Gouvernement du Pakistan. Il a recommandé que Il’Unesco 
au cas ou elle y serait conviée, accorde avx organisateurs du dit colloque 
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une assistance financiére leur permettant d’inviter des savants étrangers a 
participer aux débats. Le Comité consultatif pourrait fort bien faire coincider 
avec ce colloque sa deuxiéme session de 1957. 


Recommandation 9 


23. Le Comité a recommandé que |’Unesco s’entende avec un Etat mem- 
bre en vue d’organiser un colloque en 1958, Il a proposé que ce colloque soit 
consacré a |’étude de la salinité des sols et des eaux, et notamment 4 la puri- 
fication des eaux salines, ainsi qu’a l’utilisation des eaux salines par les plantes, 
les animaux et l’-homme. 


H — Projets spéciaux 


Recommandation 10 
24. Le Comité a pris note avec un grand intérét du voeu émis par le groupe 


de travail n° 8 (purification des eaux salines) de l’Organisation européenne 


de coopération économique dont J’objectif est maintenant atteint, tendant 
a ce qu’une organisation se charge d’assurer la liaison entre les membres 
du groupe et les autres intéressés, de diffuser les renseignements qui offrent 
de l’intérét pour le groupe, et d’organiser de temps a autre des échanges de 
vues sur les problémes abordés et les résultats obtenus. Le Comité a estimé 
que le programme de la zone aride permet de donner satisfaction a ce vceu, 
notamment par le moyen du bulletin, des rapports — inventaires et des col- 
loques. Il a recommandé en conséquence que le Directeur général accueille 
favorablement la demande du Secrétaire général de 1?OECE. II a, en outre, 
recommandé qu’un colloque sur les problemes relatifs 4 la salinité soit orga- 
nisé dans les plus brefs délais (voir recommandation 8). Ce colloque aurait 
pour objet d’adopter des recommandations touchant les travaux que le Comité 
devra entreprendre dans ce domaine. Une étroite coopération devrait étre 
assurée avec le « Saline Water Conversion Programme » réalisé aux Etats- 
Unis d’Amérique. 


J] — Publications 


Recommandation 11 


25. Considérant que le projet majeur intéressera au premier chef les pays 
du Moyen-Orient, considérant d’autre part l’intérét qu’il y a a porter, dans 
toute la mesure du possible, les renseignements concernant la zone aride 
a la connaissance de tous ceux qui pourraient en profiter, le Comité a recom- 
mandé que des mesures soient adoptées en vue de favoriser la publication dans 
les langues de la région des principaux résultats obtenus. 


REUNIONS SCIENTIFIQUES 


Recommandation 12 


26. Le Comité a entendu avec un vif intérét les exposés de ceux de ses 
membres qui ont participé au Congrés international de géographie, au Congres 
international de pédologie et au Stage d’études international sur la géographie 
réuni par l’Université musulmane d’Aligarh (Inde). Il a recommandé que ses 
membres rédigent a l’intention du Secrétariat de brefs rapports sur ces réu- 
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nions, aux fins de reproduction et de diffusion. Il serait bon que les membres 
du Comité prennent désormais l’habitude d’établir ainsi, dans la mesure du 
possible, sur les réunions scientifiques auxquelles ils participent, des rapports 
présentés selon un modéle uniforme. 


s P 
PREPARATION D’UNE CARTE AVEC LE CONCOURS DU COMITE 
DE RECHERCHES SUR LA ZONE TROPICALE HUMIDE 


Recommandation 13 


27. Le Comité a pris acte de la recommandation faite par la Réunion 
préliminaire de spécialistes des recherches sur la zone tropicale humide, tenue 
a Kandy (Ceylan), en mars 1955, touchant l’organisation d’une réunion mixte 
de représentants des Comités de recherches sur la zone aride et sur la zone 
tropicale humide en vue de préparer une carte a échelle réduite délimitant 

.dans leurs grandes lignes les régions du monde ou se posent des problémes 
intéressant les deux programmes. I] a également entendu un rapport de M. G. 
Aubert et de M. S. N. Naqvi (participants a la réunion de Kandy) sur les rai- 
sons pour lesquelles cette recommandation a été adoptée. Le Comité a invité 
le Directeur général a signaler a l’attention du Comité consultatif de recherches 
sur la zone tropicale humide les cartes homoclimatiques publiées en 1951 
dans le cadre du programme de la zone aride. Ces cartes délimitent les zones 
qui offrent un intérét particulier pour les recherches sur la zone aride. Au cas 
ou le Comité consultatif de recherches sur la zone tropicale humide jugerait 
opportun de dresser des cartes spécialement concgues a son intention et sou- 
haiterait obtenir la coopération du Comité de recherches sur la zone aride, 
celui-ci étudiera bien volontiers avec lui les problémes qui offrent un intérét 
commun. 


ENERGIE SOLAIRE 


Recommandation 14 


28. Le Comité a pris acte avec intérét : 

a) de la recommandation par laquelle le Comité consultatif international 
de la recherche dans le programme des sciences exactes et naturelles invite 
l’Unesco a maintenir la contribution qu’elle apporte au perfectionnement 
et a lVuniformisation des instruments de mesure de la radiation solaire; et 

b) de la création par l’Association internationale de météorologie d’une 
sous-commission chargée de coordonner les recherches relatives A la radia-- 
diation solaire et aux utilisations pratiques de l’énergie solaire. 

Il a exprimé le désir d’étre tenu au courant de toutes les activités entre- 
prises dans ce domaine tant par le Directeur général que par la dite sous- 
commission. 


RAPPORTS DEFINITIFS SUR LES PROJETS DE RECHERCHES 
BENEFICIANT D’UNE ASSISTANCE 
Recommandation 15 


29a Comité a pris acte de la réception des rapports définitifs concer- 
nant les projets suivants : 


a) Recherches sur la rosée, par le Professeur B. Hellstrom. 
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b) Recherches sur la rosée, par Mile I. Arvidsson. 

c) Etude phytosociologique de la flore du désert du Rajasthan, par le 
Professeur F. R. Bharucha. 

d) Recherches écologiques et éthologiques sur le chameau dans le nord- 
ouest du Sahara, par M™e G. Hautier-Pilters. 

En ce qui concerne le rapport d), le Comité estime qu’il s’agit 1a d’un 
travail de grande valeur. I invite le Secrétariat 4 s’informer s’il serait utile 
d’apporter une contribution financiére a la publication des conclusions ou a 
Vexécution des recherches complémentaires mentionnées a la fin du dit rapport. 


RAPPORTS INTERIMAIRES 


Recommandation 16 


30. Le Comité a pris connaissance avec grand intérét du rapport de M. 
P. Meigs sur l’étude des possibilités géographiques des déserts cétiers. Il aime- 
' rait recevoir le rapport définitif concernant ce projet avant sa prochaine session. 


31. Le Comité a pris connaissance avec un vif intérét du rapport du Pro- 
fesseur O. Stocker sur les recherches écologiques qu’il a effectuées en Mauri- 
tanie. 


32. Le Comité a également recu le rapport sur 1|’étude de la flore du 
Rajpoutana établie par le Professeur S. Sarup. 


33. Le Comité a écouté avec un grand intérét le compte rendu du voyage 
fait par le Professeur Emberger au Moyen-Orient au printemps de 1956. 
Il sera également heureux de prendre connaissance du rapport écrit sur le 
voyage, que le Secrétariat a été prié de reproduire et de distribuer. 


DEMANDES D’ASSISTANCE 


Recommandation 17 


34. Le Comité a examiné la requéte du Professeur Franz, de l’Université 
de Vienne. Il ne s’est pas trouvé en mesure de recommander qu’une subvention 
soit accordée sous la forme requise au’titre du voyage en Afrique Equatoriale 
francaise projeté par le Professeur Franz: aucune assurance n’est en effet 
fournie en ce qui concerne la nature de l’expédition en question, ses objectifs 
et les facilités qu’elle offrira pour la conduite de travaux scientifiques. Le 
Comité a accepté avec gratitude l’offre de M. Aubert qui s’est déclaré pret 
a aider le Professeur Franz a organiser ses travaux sur le terrain. Il a recom- 
mandé qu’une aide soit accordée au Professeur Franz quand le Secrétariat — 
aura recu de M. Aubert les plans complets du projet. 


Recommandation 18 


35. Sur la base de la recommandation 13 f) faite lors de sa 10™¢ session, 
le Comité a recommandé que le Laboratoire de climatologie de Centerton, N.J. 
(Etats-Unis d’Amérique) soit chargé de mener 4a bien |’établissement de deux 
séries de cartes climatiques du Moyen-Orient, chaque série devant comprendre 
deux feuilles. Ces cartes indiqueraient la pluviosité et l’évaportanspiration 
~ potentielle. Une fois achevées, elles seraient publiées. si possible par l’Unesco, 
en méme temps qu’une brochure expliquant quels principes ont présidé a leur 
élaboration et comment on doit les utiliser. La dépense envisagée par le Comité 
au titre de ce projet s’éléve au total a 10.000 dollars. 
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Recommandation 19 


36. Le Comité a pris acte avec un grand intérét du rapport présenté par 
l'Institut francais d’Afrique noire sur les travaux scientifiques effectués pen- 
dant l’année écoulée en ce qui concerne les enclos aménagés a Adar (Afrique 
Occidentale francaise) en vue de l’étude de la dynamique des groupeme its 
végétaux. I] n’a pas été en mesure de recommander qu’une aide soit accordée 
en vue de l’achat de l’équipement scientifique nécessaire a cette étude, son 
réglement interdisant l’emploi des fonds de l’Unesco pour Vachat d’équi- 
pement permanent, mais il a recommandé que le Directeur général alloue 
l’équivalent de 700 dollars a I’Institut pour l’aider a entretenir et a protéger 
les enclos. 


ACTIVITES PROPOSEES PAR LE POSTE DE COOPERATION 
SCIENTIFIQUE DU MOYEN-ORIENT POUR 1957 


Recommandation 20 


37. Le Comité a examiné les plans établis par le PCSMO en vue de 
participer au Troisiéme Congrés scientifique arabe qui se tiendra sans doute 
4 Damas en aoit 1957, et de fournir une assistance au Colloque sur |’ento- 
mologie du Moyen-Orient qui doit se réunir a l'occasion du cinquantiéme 
anniversaire de la Société égyptienne d’entomologie. 


38. [1 apparait que l’organisation d’une réunion consacrée spécialement 
aux problémes de la zone aride, ainsi que celle d’une exposition présentée 
dans le cadre du Congrés, seraient extrémement utiles. Le Comité pourrait 
étre représentée a cette réunion par le Professeur Thacker, qui doit se trouver 
au Moyen-Orient a l’époque indiquée, et peut-étre par un autre de ses mem- 
bres. Cette question sera examinée de nouveau au cours de la prochaine ses- 
sion; a ce moment il se peut que le Poste de coopération scientifique du 
Moyen-Orient soit a méme de communiquer des renseignements plus détaillés 
sur le Congres et sur l’organisation du Colloque entomologique du Caire. 


ASSISTANCE TECHNIQUE 


Recommandation 21 


39. Le Comité a pris connaissance avec un grand intérét du rapport 
final de la mission d’assistance technique chargée de s’occuper de la création 
d’un Institut hydrogéologique en Turquie. Il a demandé au Secrétariat de 


lui fournir, dans la mesure du possible, pour examen, des exemplaires de tous 
les rapports de ce genre. 


COOPERATION AVEC LA FAO EN MATIERE DE LUTTE CONTRE 
LES ACRIDIENS DU DESERT 


Recommandation 22 


40. Le Comité a pris connaissance avec intérét et satisfaction de la recom- 
mandation faite 4 Londres, en avril 1956, par le groupe d’experts de la FAO 
sur les mesures a long terme de lutte contre les acridiens du désert: ces experts 
ont recommandé qu’une étude écologique générale des principales aires 
de reproduction des acridiens du désert soit menée, et que la FAO invite le 
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Comité consultatif 4 y participer. Le Comité aimerait recevoir de la FAO 
des indications plus concrétes et plus détaillées sur le concours qu’il pourrait 
fournir, afin de les examiner lors de sa prochaine session. I] estime qu’il y 
aurait grand intérét a rattacher au projet majeur une entreprise collective de 
ce genre, a laquelle l?OMS pourrait elle aussi coopérer. 


CONCLUSIONS GENERALES DU COLLOQUE DE CANBERRA SUR 
LA CLIMATOLOGIE ET PLUS SPECIALEMENT LA MICRO- 
CLIMATOLOGIE DE LA ZONE ARIDE 


Recommandation 23 


41. Le Comité a constaté avec satisfaction que les discussions qui ont 
eu lieu au cours du Colloque ont mis en lumieére la nécessité d’entreprendre 
des recherches sur différents problémes concernant la zone aride, tels que 
les suivants: 

Rapports entre la composition des collectivités végétales et le régime 
hydrique et thermique; présence et utilisation des eaux souterraines, de la 
rosée et de la vapeur d’eau; maxima et minima de température et de 
luminosité qui influent sur la germination des graines; croissance et survivance 
des plantes dans les régions arides; influence exercée a cet égard par la nutri- 
tion et par la nature du sol; possibilité d’utiliser des données recueillies par 
un vaste réseau d’observateurs en ce qui concerne l’application d’indices 
climatiques; étude des populations d’insectes, d’oiseaux et d’autres animaux 
de la zone aride et de la physiologie de ces espéces; rdle et propriétés de l’iso- 
lement et de la pigmentation des couches superficielles; « proprioclimat » 
des individus; établissement des plans de laboratoires d’essai fixes et mobiles; 
distribution et effets biologiques des rayons ultra-violets. 


42. Comme I’étude de l’évaporation dépend essentiellement de la mesure 
de la perte d’eau, il est évident que pour faciliter les recherches concernant 
les divers problémes ci-dessus, il faudrait que cette mesure soit réguli¢rement 
assurée dans le plus grand nombre possible de stations climatologiques; il fau- 
drait aussi mesurer la température superficielle des lacs, des étangs et des 
cuvettes, ainsi que la vitesse du vent a une altitude suffisante pour que |’in- 
fluence de la nature du terrain ne se fasse plus sentir. Il serait préférable que 
les observations soient résumées sous la forme de totaux portant sur cinq 
ou dix jours; en ce qui concerne les radiations, la création de réseaux de postes 
d’observation utilisant des instruments robustes et d’un prix relativement 
modique est extrémement souhaitable. 


ANNEXE 
LISTE DES PARTICIPANTS 


Membres du Comité consultatif 


M. G. AUBERT (France) 

le Professeur G. V. BoGomo.Lov (Union des Républiques soviétiques socialistes) 
M. B. T. DiFrcson (Australie) 

M. M.H. Ganzi (suppléant du Professeur A. BEHNIA (Iran)) 

M. H. GREENE (Royaume-Uni) 

M. S. N. Naagvi (Pakistan) 

le Professeur R. E. G. PICHI-SERMOLLI (Italie) 


le Professeur H. O’ REILLY STERNBERG (Breésil) 
le Professeur M.S. THACKER (Inde) i 
M. C. W. THORNTHWAITE (Etats-Unis d’Amérique) 


Représentants des Institutions spécialisées des Nations Unies 


Organisation des Nations Unies pour V’alimentation et l’agriculture : M. 
R. O. WHYTE (Rome) et M. H. C. FORESTER (CSIRO, Australie) 

Organisation mondiale de la santé: M. R.N. CLarK (Genéve) 

Organisation météorologique mondiale: M. H. T. ASHTON et M. A. GIBBS 
(Bureau of Meteorology, Brisbane) 


Représentants d organisations internationales non gouvernementales scienti- 
fiques et techniques 


Conseil international des unions scientifiques: le Professeur E.S. HILLs, 
de l’Université de Melbourne 

Union des associations techniques internationales: le Professeur C. W. N. 
SEXTON, de |’Université de Melbourne. 

Union géographique internationale: le Professeur H. O’REILLY STERNBERG 

Association internationale d’hydrologie: le Professeur SEXTON 

Conseil scientifique pour l’Afrique au Sud du Sahara: M. G. AUBERT. 


Expédition Glaciologique Internationale au Groenland 


REUNION DE PARIS (14-16 décembre 1956) 


Nous ne sommes pas encore en possession du compte rendu officiel 
de cette réunion. 

Toutefois, nous considérons devoir présenter dés maintenant le texte 
suivant présenté a cette réunion par le Secrétaire de l’Association, texte qui 
précise la nature des relations existant entre l’AIHS et l?EGIG. 


Patronage de lA. 1. H.S. 

1) En accordant son patronage a VE. G.I. C., A. I. H.S. marque tout 
Vintérét qu’elle porte a l’entreprise considérée. 

L’A. I. H. S. et sa Commission des Neiges et des Glaces étant les hautes 
autorités scientifiques internationales dans le domaine de la glaciologie, 
PE. G.I. G. peut profiter de ce patronage, c’est-a-dire de la marque d’intérét 
qui lui est manifestée pour s’en prévaloir auprés des organismes capables 
de soutenir financiérement ou moralement l’entreprise ou de l’aider de facon 
quelconque. 


2) Il est nettement entendu que |’A. I. H. S. n’interviendra pas financié- 
rement dans l’entreprise. 

3) Toutefois, comme marque spéciale de l’intérét qu’elle porte A ces 
travaux, VA. 1. H.S. fera paraitre dans ses comptes rendus ou son bulletin, 
les résultats spécifiquement hydro-glaciologiques de l’expédition. 

4) VA. I. H.S. en se portant en fait garante de la valeur scientifique 
du but poursuivi par 1’E.G.1.G. et les membres de son comité directeur, 
peut en retour exiger que ces membres soient agréés par elle et par suite par 
les Comités nationaux qui, d’aprés les statuts de l’Union et de l’Association, 
choisissent les membres qui représentent chacun des pays participants. 

5) L E.G.1I.G. peut devenir une sous-commissicn de la Commission 
des Neiges et des Glaces de 1’A. I. H. S. par un vote de l’Assemblée Générale. 


18 


”— 


AVIS DE SOCIETES 


SOCIETE HYDROTECHNIQUE DE FRANCE 
199, rue de Grenelle, Paris (VIII¢) 
Invalides 13-37 


ves JOURNEES DE L’HYDRAULIQUE 
AIX-EN-PROVENCE 23/28 JUIN 1958 


La «Société Hydrotechnique de France» organise les Cinquiémes Journées 
de l Hydraulique — Il’année prochaine — du 23 au 28 juin 1958, sur le theme général : 


Turbines et Pompes hydrauliques 
qui comportera : 
— un exposé introductif sur les divers types de turbines et de pompes hydrauliques, 
et sera subdivisé en 7 questions: : 


I — Evolution de la technologie des turbines et son influence sur le Génie Civil 
II — Régulateurs des turbines 

III — Rendement des turbines et problemes de cavitation 

IV. — Exploitation des turbines 

V  — Turbines - Pompes 

VI — Pompes 


VII — Turbines de forage. 


Les séances de travail auront lieu a Aix-en-Provence les 26, 27 et 28 juin 1958 
et seront précédées d’un voyage d’études facultatif permettant de visiter des ateliers 
francais de Construction de Turbines Hydrauliques et des chantiers et aménagements 
hydroélectriques. 

Les personnes désireuses de présenter des rapports ou de prendre part aux «V™es 
Journées de |’Hydraulique » sont priées de le faire connaitre dés maintenant a la 
Société Hydrotechnique de France. 

Un résumé en dix lignes de chaque mémoire devra étre adressé en 3 exemplaires 
4 la Société Hydrotechnique de France avant le 31 octobre 1957. Les rapports retenus 
devront étre remis in extenso au plus tard le 28 février 1958. 

Les rapports seront présentés, discutés et, éventuellement, publiés en frangais. 


(Preliminary Program) 


FIRST UNITED STATES INTERSOCIETY CONFERENCE 
on 
IRRIGATION AND DRAINAGE 
Conference Theme Can man develop a permanent irrigation agriculture ? 
April 29 and 30, 1957 - Sheraton-Palace Hotel, San Francisco, Calif., U.S.A. 
Sponsored by: 


American Society of Civil Engineers (Irrigation and Drainage Division) 
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American Society of Agricultural Engineers (Soil and Water Division) 
Soil Science Society of America (Conservation, Irrigation, Drainage and Tillage 
Division) 
In cooperation with 
U.S. National Committee of the International Commission on Irrigation and Drainage 


CAN MAN DEVELOP A PERMANENT IRRIGATION AGRICULTURE? 


Monday, April 29, 9:30 AM 
Purpose of Conference. 
Teamwork in the Solution of Water Problems. 
Water in a permanent Irrigation agriculture 
Influence of Climate on Irrigation Agriculture. 
Watershed Management in Relation to Water Yield 
I. Water Yield Control Through Forest Management in Snow Pack Watersheds. 
ll. Water Yields as Influenced by Vegetation Management on Watersheds 
below Snowline. 
Ill. Importance of Phreatophytes in Water Supply. 


Monday, April 29, 2:00 PM 
Water in a permanent irrigation agriculture (Continued) 
Factors Affecting the Useful Life of Reservoirs. 
Continued Utilization of Ground Water Storage Basins. 
Can Evaporation Losses be Reduced? 
Problems of Maintaining a Satisfactory Water Quality. 
Crops in a permanent irrigated Agriculture 
Selection of Crops for Water Deficient Areas 
Cropping Practices for Maximizing Water Use Efficiency. 
Plant Diseases, Insects, and Weeds as Affected by Irrigation. 


Tuesday, April 30, 9:30 AM 
Soils in a permanent irrigation agriculture 
Soil Classification in Relation to the Irrigability of Land Areas 
Effect of Irrigation on Soil Fertility and Profile Development. 
Maintenance of Favourable Soil Structure and Adequate Water Infiltration Rates. 
Salinity and Alkalinity Hazards. 


Hydraulics of Irrigation Water Application as Affected by Soil and Other Factros. 
Drainage in Relation to a Permanent Irrigation Agriculture. 


Tuesday, April 30, —:45 PM 
Man in permanent irrigation agriculture 

Conflicting Demands for Water—Some Approaches to the Problem. 

Irrigation Organizations and Water Law As Factors in a Permanent Agriculture. 

Interstate and International Water Compacts. 

Man’s Changing Food Habits and Their Impact on Irrigation Agriculture. 

Health Problems as Affected by Irrigation Agriculture. 

Man’s Opportunities to Assure the Permanence of Irrigation Agriculture— 

A Summary. 

The programme is being organized by Mr. Robert M. Hagan, Chairman, Depart- 
ment of Irrigation, University of California, Davis, California. Member: ASAE, 
SSSA, ICID. The best available authorities in the U.S.A. in each field are being 
requested to participate in the programme. 
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INTERNATIONAL COMMISSION ON IRRIGATION AND DRAINAGE 
THIRD CONGRESS 
San Francisco, April 29 - May 4, 1957 


Under the sponsoership of the United States National Committee 
Program 
Monday, April 29: Registration 
Tuesday, April 30: Registration 
Wednesday, May 1: 10:30 a.m. Inaugural session; 2:00 p.m Technical session, 
Question No. 8. 
Thursday, May 2: 9:30 a.m. Technical session, Question No. 7. 
Friday, May 3: 9:30 a.m. Technical session, Question No. 9. 
Saturday, May 4: 9:30 a.m. Technical session, Question No. 10. 
Sunday, May 5: Local Visits 
Monday, May 6 to Friday, May 17: Study tours. 
Saturday, May 18: 11:00 a.m. Closing session at Los Angeles. 

Registration, inaugural, and technical sessions from April 29 to May 4 will 
be held at the Sheraton-Palace Hotel, San Francisco. The closing session on May 18 
will be held at Los Angeles. 

The Eighth Meeting of the International Executive Council will be held at the 
Sheraton-Palace Hotel on Tuesday, April 30, at 10:00 a.m. 

Local visits will be arranged both in San Francisco and in Los Angeles, and there 
will be a separate program for the ladies. 

The Irrigation and Drainage Division of the American Society of Civil Engineers 
will hold a Conference on April 29 and 30, and a joint meeting with the Pacific Coast 
Section of the American Society of Agricultural Engineers and the Soil Scientists 
Association of America on April 30. Delegates to the Third Congress are invited . 
to the Conference and the joint meeting. 


Questions to be discussed at the Third Congress 

As previously announced in Bulletin No. 1, the following four Questions will 
be discussed at the Third Congress: 

Question No. 7—Canal lining 

Scope of the subject—Purpose and practical aspects, particularly materials, 
technical and economic aspects, maintenance and related problems. 
Question No. 8—Soil water relationship in irrigation 

Scope of the subject—Movement of moisture in irrigated soils, special methods 
for the preservation of soil structure and fertility and the effects of soil fertility on 
duty of irrigation water. Economical use of irrigation water. 
Question No. 9—Hydraulic structures on irrigation and drainage systems 

Scope of the subject—Particularly those relating to regulation and safety, distri- 
bution and measurement of water, prevention of silting (exclusive of canals, conduits 
and river works). 
Question No. 10—Jnterrelation between irrigation and drainage 


Participants in the Congress 

‘Delegates from participating countries or from other organizations invited by 
the Central Office of the International Commission on Irrigation and Drainage, 
accompanied by their families, are invited to take part in the Congress. 

Engineers, scientists, and technical specialists from non-participating countries, 
accompanied by their families, may also take part in the Congress as individual 
participants, provided their applications are received by the central Office in New 
Delhi and by the U. S. National Committee of the Commission, not later than March I, 
1957. Applications from Western Hemisphere countries should be mailed to the 
U.S. National Committee. 
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INTERNATIONAL ASSOCIATION FOR HYDRAULIC RESEARCH 
ASSOCIATION INTERNATIONALE DE RECHERCHES HYDRAULIQUES 


SEPTIEME CONGRES - LISBONNE 
Juillet 1957 
Sujets a traiter 

Le programme définitif des sujets a traiter a été fixé comme suit: 

A — Effet d’échelle — Principes généraux, études fondamentales, applications 
aux travaux hydrauliques et aux machines. 

B — Cavitation — Principes généraux, études fondamentales, applications aux 
travaux hydrauliques et aux machines. 

C — Aydraulique des travaux de prise d’eau, de tunnels et canaux d’amenée, 
de tunnels de dérivation provisoire. 

D — Sujets libres — En ce qui concerne les sujets libres nous voulons préciser 
qu’ils ne pourront étre présentés qu’aprés approbation par le Conseil de 1’ Associa- 
tion, ceci afin d’éviter d’alourdir les travaux du Congrés et de recevoir des rapports 
dont les sujets s’écarteraient trop des themes choisis. 


COURS DE RECHERCHES SCIENTIFIQUES EN HAUTE MONTAGNE 
(Gletscherkurs) du Prof. R. Finsterwalder 


Suivant une tradition bien établie, le Professeur Finsterwalder organise un cours 
de recherches scientifiques en haute montagne. 

M. le Professeur Finsterwalder m’a communiqué, a ce sujet, lavis suivant: 
Kurs fur Hochgebirgsforschung 1957 

Der Kurs findet vom 18 bis 25 August 1957 in der Alpinen Forschungsstelle 
der Universitat Innsbruck in Obergurg!, Otztal, Tirol, statt. Er wird durch Herrn 
Prof. Dr. H. Kinzl, Innsbruck, und den Unterzeichneten geleitet. 

Das Forschungs- und Unterrichtungsprogramm des Kurses umfasst : 

Photogrametrische Aufnahme des derzeitigen Gletscherstandes in den Zungen- 
gebieten des Gurgler- und Rotmossferners sowie einzelne Aufnahmen in den Firn- 
gebieten. Gletschergeschwindigkeitsmessungen nach der photogrammetrischen Metho- 
de und mit dem Theodolit. 

Morphologische, geologische und pflanzengeographische Forschungen in den 
Gletschervorfeldern, insbesondere in den durch den Gletscherriickgang seit kurzem 
eisfrei gewordenen Gebieten, dazu geographische Aufnahmen und Untersuchungen 
in den Mordanengebieten. Studium und Aufnahme der blockgletscher. Nach M6ég- 
lichkeit geophysikalische und meteorologische Forschungen. 

Theoretische und praktische Einfiihrung in die modernen Hilfsmittel der Glet- 
scherforschung, insbesondere die photogrammetrische Aufnahme-methode, sowie 
Einfiihrung in die topographische und kartographische Aufnahme bei Forschungs- 
reisen. 

Die Zahl der Teilnehmer ist beschrankt. Alpine Erfahrung und kOrperliche 
Eignung fiir die teilweise anstrengenden Arbeiten im Gelidnde sind Voraussetzung. 
Erwiinscht sind vor allem Meldungen von solchen Teilnehmern, die sich weiterhin 
den Aufgaben der Gletscher- und Hochgebirgsforschung im In- und Ausland sowie 
der Polarforschung aktiv widmen wollen. Einige Vertrautheit mit geoditischen Mess- 
instrumenten und ihren Gebrauch ist notwendig. 

Die Kursgebiihr betragt DM 50,—; sie kann auf begriindeten Antrag erlassen 
werden. 

Le cours aura lieu du 18 au 25 aot 1957 a Obergurgl-Oetztal, Autriche. 

Pour tout renseignement s’adresser A Monsieur le Professeur R. Finsterwalder, 
Institut fiir Photogrammetrie, Arcisstrasse 21, Miinchen (Allemagne) 
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DARCY’S LAW AND THE FIELD EQUATIONS 
OF THE FLOW OF UNDERGROUND FLUIDS 


BY 


M. KING HUBBERT 


ABSTRACT 


In 1856 Henry Darcy described in an appendix to his book, Les Fontaines 
Publiques de la Ville de Dijon, a series of experiments on the downward flow 
of water through filter sands, whereby it was established that the rate of flow 
is given by the equation: 

q =— K(h.—h))/I, 


in which qg is the volume of water crossing unit area in unit time, / is the 
thickness of the sand, A, and h, the heights above a reference level of the water 
in manometers terminated above and-below the sand, respectively, and K 
a factor of proportionality. 

This relationship, appropriately, soon became known as Darcy’s law. 
Subsequently many separate attempts have been made to give Darcy’s empi- 
rical expression a more general physical formulation, with the result that 
sO many mutually inconsistent expressions of what is purported to be Darcy’s 
law have appeared in published literature that sight has often been lost of 
Darcy’s own work and of its significance. 

In the present paper, therefore, it shall be our purpose to reinform 
ourselves upon what Darcy himself did, and then to determine the meaning 
of his results when expressed explicitly in terms of the pertinent physical 
variables involved. This will be done first by the empirical method used 
by Darcy himself, and then by direct derivation from the Navier-Stokes 
equation of motion of viscous fluids. We find in this manner that: 


q = (Na*)(e/)lg —C/9) grad p] = oE , 

is a physical expression for Darcy’s law, which is valid for liquids generally 
and for gases at pressures higher than about 20 atmospheres. Here N is 
a shape factor and d a characteristic length of the pore structure of the solid, 
p and y are the density and viscosity of the fluid, o = (Nd*)(9/w) is the volume 
conductivity of the system, and E = [g —(1/¢) grad p] is the impelling force 
per unit mass acting upon the fluid. It is found also that Darcy’s law is 
valid only for flow velocities such that the inertial forces are negligible as 
compared with those arising from viscosity. 

In general three-dimensional space there exist two superposed physical 


fields : a field of force of caracteristic vector E, and a field of flow of vector q. 
The force field is more general than the flow field since it has values in all 
space capable of being occupied by the fluid. 


So long as the fluid density is constant or is a function of the pressure only, 


curl E =0, FE = —grad ® 


= er + [P. 


where 
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The field of flow, independently of the force field, must satisfy the con- 
servation of mass, leading to the equation of continuity 


div og =—fdp/dt 
where fis the porosity, and ¢ is the time. For steady motion d0/d¢ = 0, and 


div og =0. 
If the fluid is also of constant density, 
divq —0. 


The two fields are linked together by Darcy’s law, 
q=oE, 
which is physically analogous to Ohm’s law in electricity. 
Then, when Be grad MO, | 
q =o prad D. 


By means of the foregoing equations the flow of both homogeneous and 
heterogeneous fluids through porous solids becomes amenable to the same 
kind of analytical treatment as is already familiar in electrical and thermal 
conduction. 

The relation of Darcy’s work to the development of a valid theory of 
the flow of fluids through porous solids, is somewhat analogous to that of 
Faraday to the Maxwellian equations of electromagnetism. It forms a solid 
experimental foundation for such a field theory, and the errors attributed by 
various recent authors to Darcy appear upon closer inspection to have been 
those committed by the authors themselves. 


* 
* * 


In Paris in the year 1856 there was published by Victor Dalmont as a 
part of the Libraire des Corps Impériaux des Ponts et Chaussées et des Mines 
a monograph by the French engineer Henry Darcy (*), Inspector General 
of Bridges and Roads, bearing the title: 

« LES FONTAINES PUBLIQUES 
DE LA VILLE DE DIJON 
Exposition et Application 
DES PRINCIPES A SUIVRE ET DES FORMULES A EMPLOYER 
Dans les Questions 
de 
DISTRIBUTION D’EAU 
Ouvrage terminé 
par un Appendice Relatif aux Fournitures d’Eau de Plusieurs Villes 
AU FILTRAGE DES EAUX 
et 
A la Fabrication des Tuyaux de Fonte, de Plomb, de Tole et de Bitume » 


For several years previously M. Henry Darcy had been engaged in 
modernizing and enlarging the public water works of the town of Dijon, 
and this treatise, comprising a 647-page volume of text and an accompanying 
Atlas of illustrations, constitutes an engineering report on that enterprise. 

The item of present interest represents only a detail of the general work 
and appears in an appendix on pages 590 to 594 under the heading « Deter- 
mination of the Law of Flow of Water Through Sand », and pertains to a 
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problem encountered by Darcy in designing a suitable filter for the system. 
Darcy needed to know how large a filter would be required for a given quantyti 
of water per day and, unable to find the desired information in the published 
literature, he proceeded to obtain it experimentally. 

A drawing of the apparatus used is given as Figure 3 in the At/as and is 
here reproduced in facsimile as Figure 1. This consisted of a vertical iron 
pipe, 0.35 meters in diameter and 3.50 meters in length (the figure shows 
3.50 meters but the text says 2.50), flanged at both ends. At a height of 
0.20 meters above the base of the column, there was placed a horizontal 
screen supported by an iron grillwork upon which rested a column, a meter 
or so in length, of loose sand. Water could be admitted into the system by 
means of a pipe, tapped into the column near its top, from the building water 
supply, and could be discharged through a faucet from the open chamber 
near its bottom. The-faucet discharged into a measuring tank | meter square 


Appareil destine aj determiner la loi 
a travers le sable. 


de l’ecoulement de| l'eau 


Manometre 


a mercure 


Manometre 


a mercure 


Riga Fascimile of Darcy’s illustration of his experimental apparatus. (From 
Les Fontaines Publiques de la Ville de Dijon, Atlas, Figure 3) 
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and 0.50 meter deep, and the flow rate could be controlled by means of 
adjustable valves in both the inlet pipe and the outlet faucet. 

For measuring the pressures mercury manometers were used, one tapped 
into each of the open chambers above and below the sand column. The 
unit of pressure employed was the meter of water and all manometer readings 
were reported in meters of water measured above the bottom of the sand 
which was taken as an elevation datum. The observations of the mercury 
manometers were accordingly expressed directly in terms of the heights of 
the water columns of equivalent water manometers above a standard datum. 

The experiments comprised several series of observations made between 
October 29 and November 2, 1855, and some additional experiments made 
during February 18 to 19, 1856. For each series the system was charged 
with a different sand and completely filled with water. 

By adjustment of the inlet and outlet valves the water was ade to flow 
downward through the sand at a series of successively increasing rates. For 
each rate a reading of the manometers was taken and recorded as a pressure 
difference in meters of water above the bottom of the sand. The results of 
two of these series, using different sands, are shown graphically in Figure 2. 
In each instance it will be seen that the total rate of discharge increases linearly 
with the drop in head across the sand of the two equivalent water manometers. 
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FIRST SERIES Wi 
OCTOBER 29-NOVEMBER 2,1855 


i) 
je) 


FIFTH SERIES 
FEBRUARY 17-18, 1856 


DISCHARGE (LITERS PER MINUTE) 
° 


ce) 5 10 15 
DROP IN HEAD ACROSS SAND (METERS) 


Fic. 2 — Graphs compiled from Darcy’s tabular data on his experiments of 

October 29 to November 2, 1955, and of February 17-18, 1856, showing linear 
relation between flow rate and differences in heigts of equivalent water manome- 
ters. 


Darcy’s own summary of the results of his experiments is given in the 


following passage (p. 594): 
« Ainsi, en appelant e I’épaisseur de la couche de sable, s sa superficie, 
P la pression atmosphérique, / la hauteur de l’eau sur cette couche, on aura 


27 


P +A pour la pression a laquelle sera soumise la base supérieure; soient, 
de plus, P + A» la pression supportée par la surface inférieure, k un coefficient 
dépendant de la perméabilité de la couche, q le volume débité, on a 


s 
Quek 2 [h + e = ho] qui se réduit ag= k — (h + e) quand h,=0, ou lorsque 
e e 


la pression sous le filtre est égale a la pression atmosphérique. 
«Il est facile de déterminer la loi de décroissance de la hauteur d’eau h 
sur le filtre; en effet, soit dh la quantité dont cette hauteur s’abaisse pendant 


Eh nee ates ; : 
un temps dt, sa vitesse d’abaissement sera — de mais l’équation ci-dessus 


donne encore pour cette vitesse l’expression 


Mea asr yA ba 
Saat =e 8) 
dh =k dh k 
« On aura donc Tape nt Bs d’ou ae ee ot. 
k 
et Ld+ed=C— ot. 


[/ is the logarithm to the base e.] 


« Si la valeur 4) correspond au temps fy et / a un temps quelconque 1, 
il viendra 


k 
Ih +e) = Mo + €-) —— lt — Fo] (1) 
: : ge Go€ 5. 
« Sion remplace maintenant fA + e et hy + e par <a et yak il viendra 
k : 
lq = Iqo ei (t — fo) (2) 


et les deux équations (1) et (2) donnent, soit la loi d’abaissement de la hauteur 
sur le filtre, soit la loi de variation des volumes débités a partir du temps fo. 
« Si k et e étaient inconnus, on voit qu’il faudrait deux expériences pré- 


EE nea AEE ; k 
liminaires pour faire disparaitre de la seconde le rapport inconnu Pass 


Translating Darcy’s statements into the notation which will subsequently 
be used in the present paper, what Darcy found and stated was that, when 
water flows vertically downward through a sand, the volume of water Q 
passing through the system in unit time is given by 


hy h,—h 
Q=KA “>, orby— KA >; (1) 
and the volume crossing unit area in unit time by 
h — h —h, 
Q)A Sg are Ot by Kierra (2) 


where K is a factor of proportionality, A the area of cross section and / the 
thickness of the sand, and fA, and h, the heights above a standard reference 


elevation of water in equivalent water manometers terminated above and below 
the sand, respectively. 


Writing equation (2) in differential from gives 
q == —K (dh/dl) . (3) 
Soon after the publication of Darcy’s account of these experiments, the 
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relationship expressed by equations (1) to (3) became known, appropriately, 
as Darcy’s law. S 

It has subsequently come to be universally acknowledged that Darcy’s 
law plays the same role in the theory of the conduction of fluids through 
porous solids as Ohm’s law in the conduction of electricity, or of Fourier’s 
law in the conduction of heat. On the other hand, Darcy’s own statement 
of the law was in an empirical form which conveys no insight into the physics 
of the phenomenon. Consequently, during the succeeding century many 
separate attempts were made to give the statement of the law a more general 
and physically satisfactory form, with the result that there appeared in the 
technical literature a great variety of expressions, many mutually contra- 
dictory, but all credited directly or indirectly to Henry Darcy. 

It has accordingly become recently the fashion, when some of these ex- 
pressions have been found to be physically untenable, to attribute the error 
to Darcy himself. In fact one recent author, in discussing a supposed statement 
of Darcy’s law which is valid for horizontal flow only, has gone so far as to 
explain that Darcy was led to the commission of this error by restricting his 
experiments to flow in a horizontal direction. 

On this centennial occasion of Darcy’s original publication, it would 
appear to be fitting, therefore, instead of merely paying our respects to Darcy 
in the form of an empty homage, that we first establish unequivocally what 
Darcy himself did and said with respect to the relationship which bears his 
name; second, try to ascertain the generality and physical content of the 
relationship and to give it a proper physical expression; and third, attempt 
to see how this fits into a general field theory of the flow of fluids through 
porous solids in three-dimensional space. 

The first of these objectives has already been accomplished; the second and 
third will now be given our attention. 


THE PHYSICAL CONTENT OF DaRcy’s LAW 


As we have seen heretofore, what Darcy determined was that, when 
water flows vertically downward through a sand, the relation of the volume 
of water crossing unit area normal to the flow direction in unit time, to the 
thickness of the sand, and to the difference in heights of equivalent water 
manometers terminated above and below the sand, is given by the following 
equation : 

hy—h h,—h, 


Pie re ON pe Reais ; (2) 


where K is «a coefficient depending upon the permeability of the sand». 

Questions immediately arise regarding the generality of this result. 
Would it still be true if the water flowed upward through the sand? or 
horizontally? What changes would be effected in the relationship if some 
different liquid characterized by a different density and viscosity were used? 
In what manner does K depend upon the permeability of the sand, or upon 
its measurable statistical parameters such as coarseness and shape? And 
finally, what physical expression can be found which properly embodies all 
of these variables? 

The answer to most of these questions can be determined empirically 
by an extension of Darcy’s original experiment. If, for example, we construct 


an apparatus such as that shown in Figure 3, consisting of a movable cylinder 


with a rigid sand pack into which two manometers, at an axial distance / 
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REFERENCE ELEVATION 


Fic. 3 — Apparatus for verifying Darcy’s law for flow in various directions. 


apart, are connected by flexible rubber tubing, we can determine the validity 
of Darcy’s law with respect to the direction of flow. With the apparatus 
vertical and the flow downward at a total rate Q, the manometer difference 
h,—h, will have some fixed value Ah. Now, keeping Q constant and 
inverting the column so that the flow will be vertically upward, it will be 
found that Ah also remains constant. Next, setting the column horizontal, 
Ah still remains constant. In this manner we easily establish that Darcy’s 
law is invariant with respect to the direction of the flow in the earth’s gravity 
field, and that for a given Ah the flow rate Q remains constant whether the 
flow be in the direction of gravity or opposed to it, or in any other direction 
in three-dimensional space. 

This leads immediately to a generalization for flow in three-dimensional 
space. At each point in such space there must exist a particular value of 
a scalar quantity h, defined as the height above a standard elevation datum 
of the water column in a manometer terminated at the given point. The 
ensemble of such values then gives rise to a scalar field in the quantity / with 
its attendant family of surfaces, 4 = constant. In such a scalar field water 


will flow in the direction perpendicular to the surfaces, 4 = constant, and 
at a rate given by 


q = —Kgradh. (4) 
Continuing our empirical experimentation, we find that when we change 

either of the fluid properties, density or viscosity, or the geometrical properties 
of the sand, equation (4) still remains valid but the value of K changes. In 
particular, by varying one factor at a time, we find 

IO GCLD ies | 

Ko I/u, 
where ¢ is the density and wu is the viscosity of the fluid. 


Likewise, if we use a number of geometrically similar sands which 
differ only in grain size, we find that 


IKnOC ada (6) 


(5) 
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where d is a length such as the mean grain diameter, which characterizes 
the size scale of the pore structure of the sand. 

; Introducing the results of equations (5) and (6) into equation (4) then 
gives 


q — K’ dXo/y)(— grad h) , (7) 


in which K’ is a new factor of proportionality containing all other variables 
not hitherto explicitly evaluated. This remains, however, an empirical 
equation devoid of dynamical significance since there is no obvious reason 
why the flow of a viscous fluid through a porous solid should be proportional 
to a dimensionless quantity, — grad h. 

This deficiency can be eliminated when we introduce the equation relating 
the manometer height A to the dynamical quantities, gravity and pressure. 
At any point P within the flow system, characterized by elevation z and 
manometer height h, the pressure is given by the hydrostatic equation 


Pa pea 2) 5 
from which 
h =(piog) +z, (8) 
and 
— grad h = — (1/og) grad p — grad z. (9) 
Multiplying both sides of equation (9) by g then gives 
— g grad h = — (1/09) grad p — g grad z. (10) 


With the z-axis vertical and positive upward, then grad z is a unit vector 
directed upward, so that — g grad z is a vector of magnitude g directed down- 


ward. Designating this by gz, equation (10) becomes 
gs grad h —.g —(1/p) grad p, (11) 


in which each of the terms to the right represents, both in direction and 
magnitude, the force exerted upon unit mass of the fluid by gravity and by 
the gradient of the fluid pressure respectively; the fluid flows in the direction 
of, and at a rate proportional to, their resultant, —ggradh. Hence the 
dynamical factor g has evidently been concealed in the original factor K and 
must still be present in the residual factor K’. Introducing this explicitly, 
we may now write 


q = (Na?)\(p/u)(— g grad h), (12) 


in which N is a final factor of proportionality. Dimensional inspection 
shows that N is dimensionless, and a physical review indicates that no dyna- 
mical variables have been omitted, so that N must be related to the only 
remaining variable, namely the shape of the passages through which the 
flow occurs. Since shape is expressed by angular measurement and angles 
are dimensionless, [L/L], then N, must be a dimensionless shape factor whose 
value is constant for systems which are either identically, or statistically, 
similar geometrically. By identical similarity is meant similarity in the strict 
Euclidean sense : all corresponding angles equal, and all corresponding lengths 
proportional. By statistical similarity is meant that two complex geometrical 
systems which may not be identically similar on a microscopic scale are still 
indistinguishable as to shape on a macroscopic scale — for example, two 
sets of randomly packed uniform spheres. 
The term (— g grad h) can also be written in the form 


— g grad.h = — grad (gh), (13) 
where gh represents the amount of work required to lift a unit mass of water 
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from the standard datum of elevation outside the system to the height A of 
the water in the manometer. Then, since the additional work required to 
transport the water down the waterfilled tube of the manometer to its terminus 
is zero, it follows that 

gh = (14) 


is a measure of the energy per unit mass, or the potential, of the water in the 
system at the point at which the manometer is termintaed. A manometer 
is thus seen to be a fluid potentiometer, the potential at every point being 
linearly related to the manometer height 4 by equation (14). Then, if we 
let E be the force per unit mass, or the intensity of the force field acting upon 
the fluid, we have 


E =—grad ® = —g gradh = g—(1/0) gradp, (15) 


Ears 


tic. 4 — Relation of force per unit mass, E, to the primary forces 2 and 


as 


—(l/p) grad p; and of the flow vector ns to Ey 


and, by substitution into equation (12), Darcy’s law may be expressed in any 
of the following equivalent forms (Fig. 4): 
qi — (Nd?)(e/u) grad D = oE : 
q = —(Na*\(p/u) ¢ grad h = o [g — (1/0) grad p], (16) 
= —(Nd*)(9/2) g grad h = (o/0)[eg — grad p], 
where 
o = [(Nd’)\(e/)] 
is the volume conductivity fo the system. In the | i 
eueimegraities y e last of equations (16), the 
[pg — grad p] = pE = H (17) 
represents the force H per unit volume. 
When we compare Darcy’s law in the form: 


q =oE =—ograd®, 
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with Ohm’s law: 
i = 6c Ee = —o-gradV, 


where i is the current density, oe the electrical conductivity, E, the electrical 
force-field intensity, and V the electrical potential, the physical as well as the 
mathematical analogy between Darcy’s law and Ohm’s law becomes imme- 
diately apparent. 


The problem of permeability. — Now that we have achieved a complete 
physical statement of Darcy’s law, it remains for us to define what shall be 
meant by the permeability of the system. 

It wiil be recalled that Darcy stated that the factor K is «a coefficient 
depending upon the permeability of the sand» Following this it has often 
been the custom, especially among ground-water hydrologists, to define the 
permeability of a system to be synonymous with K. But, as we have seen, 
the factor K is the lumped parameter, 

K = (Na’)(o/u) g , 
comprising the geometrical properties of the sand, the dynamical properties 
of the fluid, and even the acceleration of gravity. Consequently, if K is taken 
as a measure of the permeability, it will be seen that the same sand will have 
different permeabilities to different fluids. 

During recent years there has been a convergence of opinion toward 
the conclusion that permeability should be a constant of the solid indepen- 
dently of the fluids involved. If we accept this view, then it is seen that the 
only property of the solid affecting the rate of flow is the geometrical factor, 

k=Nad’, 
which we may accordingly define to be its permeability. Then, since N is 
dimensionless and d is a length, it follows that the dimensions of permeability 
are [L”]; and in any consistent system of units, the unit of permeability is the 
square of the unit of length. 

In practice the magnitude of this quantity, for a given porous solid, is 
determined hydrodynamically by flowing a liquid through the solid, and mea- 
suring all variables except k, and then solving Darcy’s law for k: 


Kee NG? = esse cee —- 
eg — grad p 
__ which, when the vector quantities are resolved into their components in the 
_ flow direction s, becomes 


qv. 
Saved ek alas e: Rr) 
It is found in this manner that for randomly packed, uniform spheres 
of diameter d, the value of the shape-factor N is approximately 6 x 10~‘, 
- Then, for a pack of uniform spheres of any size, the permeability will be 
approximately 
k=(6-<10-*).a"* 
If d, for different packs, is allowed to vary from about 10—* to 10~—' cm, 
corresponding to the approximate range of grain sizes from fine silts to coarse 
sands, the permeability will vary from about 10-1? to 10—°cm?, which is 
also approximately the range of the permeabilities of the corresponding clastic 
sediments. 
In view of the fact that magnitudes of permeabilities of rocks are remote 
from that of the square of any unit of length in commun use, there is some 
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advantage in having a practical unit such that most measured values fall 
within the range 1 — 10,000 practical units. If such a practical unit is to 
fit into a consistent system of measurement without awkward conversion 
factors, then it must also be a submultiple of the fundamental unit of the form: 
1 practical unit = 10—” fundamental units. 
In the cgs system with the fundamental unit the (centimeter)?, the optimum 
value of the exponent would be about, 12 or 


1 practical unit = 10—'* cm*. 
Regrettably the unit of permeability used almost universally in the 


petroleum industry, for which the name «darcy» has been preempted, was 
defined originally in terms both of an incomplete statement of Darcy’s law :2:3 


| ee 
i> cata rice (19) 


and an inconsistent system of measurement. The permeability k is defined 


to be 1 darcy when g = 1 (cm*/cm?)/sec, 1 = 1 centipoise, and dp/ds = 1. 


atmosphere/cm. 

In the complete Darcy’s law of equation (18), the factor pgs, except 
when the motion is horizontal, is of comparable magnitude to dp/ds, and 
so cannot be ignored. Consequently, since it is not practical to measure 
egs in atmospheres/cm, it follows that permeabilities, expressed in darcys, 
cannot be used in a proper statement of Darcy’s law without the insertion 
of a numerical factor to convert ogs from cgs units into atmospheres/cm. 

The only alternative is to convert permeabilities expressed in darcys 
into the cgs unit, the cm®. For this conversion 

1 darcy = 0.987 x 10—® cm?, 


which is within 1.3 percent of the submultiple, 10—* cm’. 

Few permeability measurements are accurate to within 1.3 percent and, 
when several specimens from the same formation are measured, the scatter 
is much greater than this. Consequently, for all ordinary computations, 
the approximate conversions : 

1 darcy = 10—' cm?, 
1 millidarcy ~ 10—11.cm?, 


are more accurate than the permeability data available; though if the data 
warrant it, the more precise conversion can of course be used. 


DERIVATION OF DARCyY’s LAW FROM NAVIER-STOKES EQUATION 


Having thus achieved the desired generalization and a proper physical 
statement of Darcy’s law by an extension of the empirical method which 
Darcey himself employed, let us now see if the same result can be derived directly 
from the fundamental equation of Navier and for the motion of a viscous fluid. 


Macroscopic and microscopic scales.—In order to do this we must first 
distinguish between the two size scales, the macroscopic and the microscopic 
on which the phenomena considered are to be viewed. : 
The macroscopic scale, which is the one we have been using thus far, 
is a scale that is large as compared with the grain or pore size of the porous 
solid. On this scale the flow of a fluid through a porous solid is seen as a 
continuous phenomenon in space. However, when we are dealing with 
macroscopic quantities which have particular values at each point in space, 
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but which may vary with position, it is necessary for us to define more clearly 
what is meant by the value of a macroscopic quantity at a given point. 

This can be illustrated with the concept of porosity. Suppose that we 
are interested in the porosity at a particular point. About this point we take 
a finite volume element AV, which is large as compared with the grain or pore 
size of the rock. Within this volume element the average porosity is defined 
to be 


oS  Aafoe (20) 


where AV; is the pore volume within AV. We then allow AV to contract 


about the point P and note the value of f as AV diminishes. If we plot if 
as a function of AV (Fig. 5), it will approach smoothly a limiting value as 
AV diminishes until AV approaches the grain or pore size of the solid. At 
this stage f will begin to very erratically and will ultimately attain the value 
of either 1 or 0, depending upon whether P falls within the void or the solid 
space. 


(9) AV —~> 


Fic. 5 — Method of defining point values of macroscopic quantities illustrated with 
the porosity /f. 


However, if we extrapolate the smooth part of the curve of +3 versus AV 
to its limit as AV tends to zero, we shall obtain an unambiguous value of f 
at the point P. We thus define the value of the porosity fat the point P to be 


extrap lim AV; 
FP)= av_so AV’ 


where extrap lim signifies the extrapolated limit as obtained in the manner 


(21) 


- just described. 


By an analogous operation the point value of any other macroscopic 
quantity may be obtained, so that hereafter, when such quantities are being 
considered, their values will be understood to be defined in the foregoing 
manner, and we may state more simply: 
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lim 


or 
lim is 
OUP ES NG aa ty se (AS), > (22) 
or 
lint 
QP) = qos DAD), 
where the quantity of interest is a function of a volume, an area, or a length, 
respectively. 


The microscopic scale, on the contrary, is a scale commensurate with the 
grain or pore size of the solid, but still large as compared with molecular 
dimensions or of the motional irregularities due to Brownian or molecular 
movements. 


Microscopic equations of motion——Let us next consider the steady, 
macroscopically rectilinear flow of an incompressible fluid through a porous 
solid which is macroscopically homogeneous and isotropic with respect to 
porosity and permeability. We shall then have the fluid flowing with a constant 


macroscopic flow rate q under a constant impelling force per unit of mass E, 
and in virtue of the isotropy of the system, we shall have 

os > : 

Ue et 28 (23) 


where o is an unknown scalar whose value we shall seek to determine. 
Then, choosing x-, y-, and z-axes, 


q = (qx + jqy tka, | 


(24) 
PAT p Eo a 
and 
qx = GEx, 
qy =OEy, > (25) 
qz =OE; . 


where /, /, and k are unit vectors parallel to the x-, y-, and z-axes, respectively, 
and the subscripts signify the corresponding scalar components of the vectors 
Further, there will be no loss of generality, and our analysis will be so- 


me what simplified, if we choose the x-axis in the macroscopic direction of 
flow. Then 


= da. Pal Bs 
and (26) 
dy = Gz = 0; Ey = BE, = 0. 

Next, consider the microscopic flow through a macroscopic volume 
element AV of sides Ax, Ay, and Az. The void space in such an element will 
be seen to be an intricately branching, threedimensional network of flow 
channels, each of continuously varying cross section. A fluid particle passing 
through such a system will follow a continuously curving tortuous path. 
Moreover, the speed of the particle will alternately increase and decrease 
as the cross section of the channel through which it flows becomes larger or 
smaller. Such a particle will accordingly be seen to be in a continuous state 
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of acceleration with the acceleration vector free to assume any possible direc- 
tion in space. 


Consider now the forces which act upon a small volume element dV 
of this fluid. By Newton’s second law of motion 
dma = XdF (27) 


* * =: 5 wal brace 
where dm is the mass of the fluid, a the acceleration, and SdF is the sum of 
all the forces acting upon the fluid contained within dV. There are many 
ways in which these forces may be resolved, but for present purposes it will 


be convenient to resolve them into a driving or impelling force dFg and a 
resistive force arising from the viscous resistance of the fluid element to 


deformation, dF». Equation (27) then becomes 
dma =dFa + dFo. (28) 
_ By the principle of D’Alembert we may also introduce a force 
~ dFy = — dma, which is the inertial reaction of the mass dm to the acceleration 
a, and with this substitution equation (28) becomes 


GPa de, dks = 0. (29) 


Of these forces, dFa, which is imposed from without and does not depend 
primarily upon the motion of the fluid, may be regarded as the independent 


variable. The forces dF, and dF, both owe their existences to the fluid motion, 
and their effect is to impede that motion. 

The relation of the separate terms of equation (29) to the externally 
applied forces and the fluid motion are given by the equatoin of Navier and 
Stokes, which in vector form may be written as follows: 


elg —(1/9) grad p]dV = 9 (Dv/Dt) dV —u [yy + 1/3) Vy -dV, 30) 


in which y is the microscopic velocity, and p = — (1/3)(om + oy + oz) is 
the microscopic pressure at a point. The o’s are normal components of 
microscopic stress. 


The expression Dy/Di is the total derivative with respect to time of the 
velocity y, and is equal to the acceleration a. This can be expanded into 


—> dv — ae 
Dy/Di = — v.vY, 
y/ ¥ ee AY 


in which the term dv/dt signifying the rate of change of the velocity at a parti- 


cular point is zero for steady motion. The expression v. v, in the last term 
of equation (30), is the divergence of the velocity; and for an incompressible 
fluid this also is zero. 

Since the flow being considered is the steady motion of an incompressible 
fluid, equation (30) simplifies to 


o[g — (1/9) grad p] dV — ep (v. yr) dV —uly*») dV. (31) 


Transformation from microscopic to macroscopic equations of motion.— 
Equation (31) expresses the relation of the fluid velocity and its derivatives 


in a small microscopic volume element to the applied force dFa acting upon 
that element. If we could integrate the three terms of this equation with 
respect to the volume, over the macroscopic volume AV, and then convert 
the results into equivalent macroscopic variables, our problem would be solved. 
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In fact the integration of the first two terms presents no difficulty. The total 
driving force on the fluid content of the volume AV, as obtained by integrating 


the microscopic forces dFqg, is: 
Fa = | dFa = 5 (eg — grad p) dV = (pg — grad p) fAV, (32) 
fAV 
where grad p is the volumetric average of the microscopic grad p over the 
fluid volume fAV. 
From the macroscopic equations, the driving force Faq is given by: 
Fa = e fAV E = ele — (1/p) grad pl fAV . (33) 
Then by combining equations (32) and (33), 


F (og —grad p)dV = (eg —grad p) fAV=¢ fAV[g—(1/9) gradp], (34) 
Vv 
from which it is seen that the macroscopic gred p is equal to the volumetric 


average, grad p, of the microscopic grad p. 
Integrating the inertial term: 


F,=JdFa=e J v.yvdv 
fAV 


=io \ (wdu/dx + vdu/dy + wdu/dz) dx dy dz 
fAV 
Es (35) 
sah © j (udv/dx + vdv/dv + w dv/dz) dx dy dz 
fAV 
ap kp 


f (udw/d0x + vdw/d¥ -+ wdw/dz) dx dy dz. 
SAV 


Since there is no net gain in velocity with macroscopic distance, each 
separate integral of the expanded form of equation (35) is equal to zero, 
and we obtain for the volume element AV, 


Pp lars oe (36) 
Then, in virtue of equation (36), 
\dFo = —w J yvdvV =— Fy (37) 
fAV : 
Ordinarily the evaluation of 
j vy dV 
f AV 


would require a detailed consideration of the geometry of the void space 
through which the flow occurs and of the flow field within that space. This 
difficulty can be circumvented, however, and the integral evaluated except 


for a dimensionless factor of proportionality, provided the flow field is kinema- 
tically similar for different rates of flow. 


Criteria of similarity. Consider two flow systems consisting of two 
geometrically similar porous solids through which two different fluids are 
flowing. The criterion of geometrical similarity is that if /, and J, are any 


38 


corresponding lengths of the two systems, then for every pair of such lengths 
13/1, =a ke =S const. (38) 


The criterion of kinematic similarity is that if y, and v, are the velocities 
at corresponding points in the two systems, the two velocities must have 
the same direction and their magnitudes the ratio 


Vo/V, = Vr = const . 


Then, since the forces dF, and dF» acting upon a fluid element are each 
determined by the velocities and the fluid density, or viscosity, and dFy is 


> > 
determined by dF, and dF,, if the fluid motions of the two systems are kine- 
matically similar, all corresponding forces will have the same directions and 
their magnitudes the same ratio (Fig. 6). 


FiG. 6 — Microscopic views of two dynamically similar flow systems 


Thus 
(dFa)2 (dF »)s (dFa)2 
(dFa), bx (dF); (dFa)1 ; 


Since only two of the three forces are independent, we need to consider the 
ratios of only the first two, and, by reciprocation, 


dFa dFa 
ALD py (ee 40 
ee aie (40) 


indicating that for each system the ratio of the inertial to the viscous force 
must be the same. 


(39) 


From equation (31) 
(dFa)2 2 (Vv. Yv)a Ve | 


+ ; (41) 
(dFa), Qi (Vv. YW dV, 
(dFv)» eS Ue (VV)2 dV (42) 


(dFy): Ui (Vv) dvi 
In equation (41) p. vy expands into the sum of a series of terms, each 
of the form uw du/dx, which is a velocity squared divided by a length. In 
equation (42) V7’v expands into a series of terms, each of the form 07u/0dx°, 
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which is a velocity divided by the square of a length. Then, since the ratios 
of all corresponding velocities and of all corresponding lengths of the two 
systems are constant, we may choose any suitable velocity and any convenient 
length. We accordingly choose for the velocity the macroscopic flow rate q 
whose dimensions are [L* L—? T—], or [LT—"]. For the characteristic length 
we choose d, which must be some convenient statistical length parameter 
of the microscopic geometry of the system. With these substitutions the 
first two terms of equation (39) become 


02 2” ds" {Lo G2 dz 


eat qd [1 qi d,’ Co 
and, by reciprocation, equation (40) becomes 

O22" d,” 01 91° dy’ 

Le Jo Ay in: {hy qd, 
or 

q2 ds -. ad, (44) 


o/P2 YrlPr 
The dimensionless quantity (qd)/(u/e) is the Reynolds number R of the 
system, which, as seen from its derivation, is a measure of the ratio of the 
inertial to the viscous forces of the system. Our criterion for kinematic 
similarity between the two systems thus reduces to the requirement that 
Ree — lite = (45) 
Now let us specialize the two systems by making 


d, =a, 061 =.01, Ue =a; 


which is equivalent to requiring the same fluid to flow through the same porous 
solid at velocities, g. and g,. However, when these values are substituted 
into equation (44), we obtain 


Gqz2=— M1, 
indicating that, in general, when the same fluid flows thourgh a given porous 
solid at two different rates, the resulting flow fields cannot be kinematically 
similar. 

However, since dFq is proportional to q? and dF, to g, then as gq is de- 
creased dF, diminishes much more rapidly than dF,. Consequently there 
must be some limiting value of g = q*, or of R = R*, at and below which 
the inertial force dFa is so much less than the viscous force dF, that the effect 


of the former is negligible as compared with the latter. For flow in this 
domain we may then write 


dod Fak 
and the force ratios become 
(dFy)» uh (dFa)> 
(dFy), i (dFa), : coe) 


whereby kinematical similarity is maintained for all rates of flow g < q*. 


Integration of viscous forces.—With this result established let us now 
return to the integration of dF» over the volume element SAV. 


) d&=—uli J yuavejs { ywavek | ywavl. (47 
fAV fAV fAV pay” 


In virtue of the fact that, by our choice of axes, dy and qz are both zero 
40 
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and there is no net flow in the y- or z-direction, the last two integrals to the 
right are both zero, and equation (47) simplifies to 


§ dFw=—ip) yuav. 48 
fAV Mav” on 
From our earlier discussion, so long as the flow remains kinematically 
similar for different rates, the quantity V’u, which is a velocity divided by 
the square of a length, is related to the macroscopic parameters by 
Vu = «a (q/d’) , (49) 
where « is a dimensionless constant of proportionality for the element dV 
but has a different value for each different element of volume. 
Substituting equation (49) into equation (48) then gives 
| dFy =—iu(qid) J adv =— "4 fr 
a eed uw (q/ ae NeOY > (50) 
where 1/N is the average value of « over fAV. 
Substituting this result into equation (37), we obtain 


—- 


elg — (/p) grad p] (AV = —J dF. ee tay , 


or 


q = Nd? (9/u) [g — (1/9) grad p), (51) 


which is the derived Darcy’s law in the same form as equation (16) deduced 
earlier from empirical data. 


Discussion of Darcy’s law.—The direct derivation of Darcy’s law from 
fundamental mechanics affords a further insight into the physics of the pheno- 
mena inyolved over what was obtainable from the earlier method of empirical 
experimentation. It has long been known empirically, for example, that 
Darcy’s law fails at sufficiently high rates of flow, or at a Reynolds number, 
based on the mean grain diameter as the characteristic length, of the order 
Ola Rj==-12= 

At the same time one of the most common statements made about Darcy’s 
law has been that it is a special case of Poiseuille’s law; and most efforts 
at its derivation have been based upon various models of capillary tubes or 
of pipes. It also has been known since the classical studies of Osborne Rey- 
nolds® in 1883 that Poiseuille’s law fails when the flow makes the transition 
from laminar to turbulent motion, so the conclusion most often reached as 
to the cause of the failure of Darcy’s law has been that the motion has become 
turbulent. 

From what we have seen, this represents a serious misinterpretation and 
lack of understanding of Darcy’s law. In the Darcy flow each particle moves 
along a continuously curvilinear path at a continuously varying speed, and 
hence with a continuously varying acceleration; in the Poiseuille flow each 
particle moves along a rectilinear path at constant velocity and zero 
acceleration. Therefore, instead of Darcy’s law being a special case of 
Poiseuille’s law, the converse is true; Poiseuille’s law is in fact a very special 
case of Darcy’s law. Another special case of Darcy’s law is. the rectilinear 
flow between: parallel plates. 

Consequently deductions concerning the Darcy-type flow made from 
the simpler Poiseuille flow are likely to be seriously misleading. The deduc- 
tion that the Reynolds number at which Darcy’s law fails is also the one at 
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which turbulence begins is a case in point. We have seen that the cause 
of the failure of Darcy’s law is the distortion that results in the flowlines 
when the velocity is great enough that the inertial force becomes significant. 
This occurs at a very slow creeping rate of flow which, for water, has the 
approximate value of 

ae (I< 10" cm seo) (ay 


Giese 


mG 


corresponding to R* = 1, when d is the mean grain diameter. 

Thus when d = 10-2.cm Darcy’s law fails at a flow rate q = | cm/sec. 

Since this represents the threshold at which the effects of inertial forces 
first become perceptible, and since turbulence is the result of inertial forces 
becoming predominant with respect to resistive forces, it would be inferred 
that the incidence of turbulence in the Darcy flow would occur at very much 
higher velocities or at very much higher Reynolds numbers than those for 
which linearity between the flow rate and the driving force ceases. That 
this is in fact the case has been verified by visual observations of the flow 
of water containing a dilute suspension of colloidal bentonite through a trans- 
parent cell containing cylindrical obstacles. This system, when observed 
in polarized light, exhibits flow birefringence which is stationary for steady 
laminar flow but highly oscillatory when the motion is turbulent. Obser- 
vations of only moderate precision indicate that the incidence of turbulence 
occurs at a Reynolds number of the order, of 600 or 700, or at a flow velocity 
of the order of several hundred times that at which Darcy’s law fails. 


Examination of the shape-factor N.—In both procedures used thus far, 
the factor N has emerged simply as a dimensionless factor of proportionality 
whose magnitude is a function of the statistical geometrical shape of the 
void space through which the flow occurs. For systems which are either 
identically or statistically similar geometrically, N has the same value. Beyond 
this we have little idea of the manner in which N is related to the shape or 
of what its numerical magnitude should be, except as may be determined 
by experiment. Let us now see if the value of N, at least to within an order 
of magnitude, can be determined theoretically. 

Since we have already seen in equation (50) that 1/N = a, where « is 
the factor of proportionality between the macroscopic quantity g/d® and the 
microscopic quantity Vj°u, it follows that in order to determine the magnitude 
of N we must first determine that of the average value of yu. For this 
purpose, with the fluid incompressible, the macroscopic flow parallel to the 
x-axis, and the inertial forces negligible, only the x-component of the Navier- 
Stokes equation (31), 


08» — dp/dx = u( 


OU oe One Oe 5 
Ix? ay? 2? iy b) (5 ) 


needs to be considered. When this is integrated with respect to the volume 
over the fluid space fAV, it becomes 


op es 0? da? 2 
(ge —dp/dx) fAV = out oa dv + aE 7 dV + : ~ dV | . (53) 

Zz 

fAV fAV fAV 

Of the three integrals to the right, the first, which represents the expansion 
of the fluid in the x-direction, is zero; and from symmetry, the last two, both 
being the integrals of derivatives with respect to axes at right angles to the 
flow, are equal to each other. In consequence, equation (53) is reduced to 
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the simpler form: 


(pg —dp/dx) fAV = —2p\(d2u/d»*) dV = —2u.Q2u/dy2) fAV. 
in which 0?u/dy? is the average value of 0*u/dy* throughout the macroscopic 
volume element. Solving this for d?u/dy? then gives 


(0*u/dy*) = —(1/2u)(pg2 — dp/dx) . (54) 

Our problem now reduces to one of attempting to determine the average 
value of 07u/dy* for the system in terms of the kinematics of the flow 
itself. If we extend any line through the system parallel to the y-axis, this 
line will pass alternately through solid and fluid spaces. At each point on 
the line in the fluid space, there will be a particular value of the x-component 
u of the velocity, which will be zero at each fluid-solid contact, but elsewhere 
will have finite, and usually positive values, giving some kind of a velocity 
profile across each fluid gap. If this profile for each gap could be determined, 
then we could also compute 0°u/dy*? at each point on the line and thereby 
determine 0°u/dy? for the line, which, in a homogeneous and isotropic system, 
would also be the average value for a volume. 

To attempt to do this in detail would be a statistical undertaking beyond 
the scope of the present paper. As a first approximation, however, we may 
simplify the problem by assuming: of 

1. That all the gaps are equal and of width 2A, where A is the average 
half-width of the actual gaps. 

2. That through each gap the velocity profile satisfies the differential 
equation 


o*u/dy? = d07u/dy? = —C. (55) 

3. That the total discharge through the averaged gaps is the same as 
that through the actual gaps. 

The velocity profile for the averaged gaps can then be obtained by inte- 
grating equation (55) with respect to y. Taking a local origin of coordinates 
at the middle of the gap, and integrating equation (55) twice with respect 
to y, we obtain 

u=—Cy?/2+Ay+B, (56) 
in which A and B are constants of integration. Then, supplying the boundary 
conditions, u = 0 when y = + i, gives 


O'= = C7 /2 + A+B 


from which 
A =0; B =Ci?/2. 
Substituting these into equation (56), we obtain 
“= {C/2)0°— y*) NSD) 


as the equation of the parabolic profile of the averaged velocity across the 
averaged gap. e 
The mean value, uw, of u across this gap is given by 


s 2 yk Cashes en 
a= (/) { udy == (Q — y*)dy = CX /3. (58) 
0 2h 


Then, replacing C by (1/2u)(eg2 — dp/dx) from equations (55) and (54), 
we obtain 


[gx — (1/9)dp/dx] . (59) 


a| 
e lo 


43 


This can be converted into terms of the macroscopic velocity, qu, by 
noting that for a macroscopic length of line / normal to the flow direction 


dal = 2u DA > 
or 
dx = 2u(XA/) = uf (60) 
where f is the porosity. 
With this substitution equation (59) becomes 


ge HLS © [en — (1/9) dp/pal (ev 


or, more generally, 


a 
g=!% . © t¢—(/p)arad al, (62) 
6 wp 

which is a statement of Darcy’s law that is valid to the extent of the validity 
of the averaging approximation used in its derivation. 

In this, it will be noted, that the geometrical factor ( f/6)22 represents 
the permeability, so that 

k (flO) = Ny”, 
or 
Nye 16:3 (63) 

where Ny, is the shape factor corresponding to ) as the characteristic length 
of the system. 


Determination of }—The mean half gap-width ) along a linear traverse 
can be determined in either of two ways. The most obvious way is by direct 
observation by means of micrometer measurements along rectilinear traverses 
across a plane section of the porous solid. 

Of greater theoretical interest, however, is an indirect method due to 
Corrsin °. Instead of a line, let a rectangular prism of cross-sectional area 
5”, where 5 can be made arbitrarily small, be passed through a porous solid 
which is macroscopically homogeneous and isotropic. This prism will pass 
alternately through solid segments and void segments. Let nm be the number 
of each which is traversed per unit length. Then the number of intersections 
with the solid surface per unit length will be 2m, and if a is the average area 
of the solid surface cut out by the prism at each intersection, the total area 
per unit length, dB, will be 

dB =2na. (64) 


If a parallel family of such prisms is made to fill all space, the number 
per unit area perpendicular to the axis of the prisms will be 1/5?, and the 
solid surface intersected per unit volume will be 

8 = 2n (a/5*). (65) 


‘ i addition, the total length of the void spaces per unit length of line 
will be 


Ore 
or 


n=fiQh): (66) 
Substituting the value of n from equation (66) into (65) and solving for ) then 
gives 
f 


a: 


i= (67) 
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Since f and 8 can be measured, the value of 2 could be determined if 
a/5® were known. 

The latter can be determined in the following manner: A homogeneous 
and isotropic distribution of the internal surface S, inside a macroscopic 
space, implies that if all equal elements dS of the surface were placed without 
rotation at the center of a reference sphere, their normals would intersect 
the sphere with a uniform surface density. Then, with the normals fixed 
in direction, if the surface elements were all moved equal radial distances 
outward, at some fixed radius they would coalesce to form the surface of a 
sphere. If the prisms of cross-sectional area 52, parallel to a given line, were 
then passed through this sphere, the average value of a/S2 would be 

it ee 
ald Sy: (68) 


and, when the summation includes one whole hemisphere, 


as? pipes area of hemisphere 2rr? 


= = 2 69 
area of diametral plane zr? be 


Introducing this result into equation (67) then gives 


A 27/6: (70) 

A method for measuring & has been described by Brooks and Purcell ’, 

but for present purposes the data on randomly packed uniform spheres, for 

which % can be computed, will suffice. For such a system, with n spheres 
per unit volume, 


area of spheres n.4r? ss Ot (1) 
~ total volume a Are’® fy Be Bak eee 
PE get ae) 


where d is the sphere diameter. 

Then, introducing this result into equation (70), we obtain for packs 
of uniform spheres 

he 


sped et! i 
a(t) 


f? 
Sep oui NSE get 
iio say) 


m (72) 


and 
(73) 


Comparison with experimental data.—In order to compare the value of 
Nx, corresponding to the characteristic length , with Na, corresponding to 
the sphere diameter d, we must first establish the relation between Nx and Ng. 
This can be done by noting that, by definition, 


Nya” = Nad? =k, 
where k is the permeability of the system. Consequently 
Nx = Na (@?/0") . (74) 


Then, introducing the value of al” from equation (73) into equation (74) 
gives for a pack of uniform spheres 


N= ler a Na. (75) 

The value of Ny for well-rounded quartz sands, screened to nearly 
uniform sizes, has been determined by the author’s former research assistant, 
Jerry Conner. Using packs of different uniform sands with mean grain 
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diameters from 1.37 X 10-? to 7.15 x 10-? cm, Conner made seven inde- 
pendent determinations of Na. The average value obtained was 
No = 610 xo10—7 (76) 
with individual values falling within the range between 5.3 x 104 and 
lest Wiles 
Conner did not determine the porosity, but the average value of the 
porosity of randomly packed, uniform spherical glass beads found by Brooks 
and Purcell? was 0.37. Inserting this value of f, and Conner’s value of Na, 
into equation (75), then gives for the equivalent experimental value of Nz: 
; Na. = 26,2 Na =i5i x 10s. (77) 


Comparing this experimental results with the approximate theoretical 
result of equation (63), it will be seen that 


Nx (theoretical) _ f/6 ve 
N>, (observed) 1.57. x 10-2 


Since our object at the outset was merely to gain some insight into the 
nature of the shape-factor N occurring in Darcy’s law, no particular concern 
is to be felt over the discrepancy in equation (78) between the observed and 
the theoretical values. All that this really indicates is that the system of avera- 
ging required should be better than the oversimplified one actually used. 
For a more complete analysis, account needs to be taken of the frequency 
distribution of the half gap-width A, and also of the functional relation between 
u and >. The fact that our approximate analysis yields a result in error by 
only a factor of 4 makes it appear promising that if account is taken of the 
variability of A and of uw as a function of A, much better approximations may 
be obtainable. 


= 4. (78) 


Darcy’s law for compressible fluids —Our analysis thus far has been rest- 
ricted to the flow of incompressible fluids for which the divergence term, 


= (1/3).2 VV. could be eliminated from the Navier-Stokes equation. 
For the flow of a compressible fluid, this term must be retained, and with 
the flow parallel to the x-axis, 


AOE OL Meee tel me Oey ow 
estes I (5 Sa Alege se tae 79 
Bes ea Pa WS ys (ae o) (79) 
of which the last two terms are equal to zero, leaving 
4 DP Man Teas OL 
Fa=—u f (22 a) 
z UM ( 30x? ze oy? a: oz? oye 


; 


or 


Fa = — 2p [(2/3)(d02u/dx2) + (02u/dy2)] fAV . (80) 


Here, OPuldx? represents the gradient of the divergence of the velocity 
in the x-direction, or the rate of the fluid expansion. Should this term be 
of the same order of magnitude as d°u/dy?, and if the flow of a gas through 
the porous system is otherwise similar to that of a liquid, then the viscous 
resistance to a gas should be greater than that for a liquid of the same 
viscosity. 

To compare the two terms d2u/dx? and 07u/dy?, it will be noted that 
each is of the form: velocity/(length)?. We have already seen that 


dufdy? = — 3 ulr? & —- gr, 
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indicating that the magnitude of this term is determined by the fact that large 
variations of u in the y-direction take place within the width of a single pore. 
Comparable variations of u in the x-direction, however, due to the expansion 
of the fluid, occur only in fairly large macroscopic distances. Consequently 
we may write 


07u/dx? & 0°g/dx? & g/P, 

where / is a macroscopic distance. The ratio of the two terms is accordingly 

07u/dx? r? 

— 2 81 

On (31) 
Then, since /?> == it follows that the additional frictional drag caused by 
the divergence term is negligible, and this term may be deleted from the 
equation. 

We conclude, therefore, from this approximate analysis, that Darcy’s 
law in its differential form is the same for a gas as for a liquid, provided that 
the flow behavior of a gas in small pore spaces, other than expansion, is similar 
to that of a liquid. 

It has’ been conclusively shown, however, by L. J. Klinkenberg * that 
the two flows are not similar, and that, in general, kg, the permeability to gas 
based on the assumed validity of Darcy’s law for gases, is not equal to ke, 


the permeability to liquids; and, in fact, is not even a constant. 
In the case of the flow of a liquid through small pores, the microscopic 


@ 
| 


velocity v becomes zero at the fluid-solid boundary; for gas flow, on the con- 
trary, there exists along the boundary a zone of slippage of thickness 5, which 
is proportional to the length of the mean-free path of the molecules. Conse- 
quently the gas velocity does not become zero at the boundaries, and the 
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Fic. 7 — Variation as a function of 1/p of the apparent permeability of a given solid 
as determined by the different gases. The value of 2.75 millidarcys, as (1/p) 
tends to 0, differs but slightly from that of 2.55 millidarcys obtained using a 
liquid (after L.J. Klinkenberg, API Drilling and Production Practice, 1941) 
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frictional resistance to the flow of gas is less than that for a liquid of the same 
viscosity and macroscopic velocity. es 

Since 5 is proportional to the mean-free path, it is also approximately 
proportional to I/p. Consequently when the gas permeability, kg, of a given 
porous solid is determined with the same gas at a number of different mean 
pressures, the resulting values of kg, when plotted as a function of 1/p, give 
a curve which is approximately linear with 1/p. Moreover, different gases, 
having different mean-free paths, give curves of different slopes. The limiting 
value of kg as 1/p—> 0, or as p becomes very large, is also equal to kj, the 
permeability obtained by means of a liquid (Fig. 7). 

In view of this fact it is clear that, in general, the flow of gases through 
porous solids is not in accordance with Darcy’s law. However, from Klin- 
kenberg’s data, at pressures greater than about 20 atmospheres (2 x 10° 
dynes/cm?, or 300 psi), the value of kg differs from ky by less than 1 percent. 
Therefore, since most oil and gas reservoir pressures are much higher than 
this, it can be assumed that gases do obey Darcy’s law under most reservoir 
conditions. 


FIELD EQUATIONS OF THE FLOW OF FLUIDS THROUGH PoROUS SOLIDS 


The establishment of Darcy’s law provides a basis upon which we may 
now consider the field equations that must be satisfied by the flow of fluids 
through porous solids in general, three-dimensional space. This problem 
is complicated, however, by the fact that the fluids considered may be either 
of constant or of variable density; that one or several different fluids, either 
intermixed or segregated into separate macroscopic spaces, may be present 
simultaneously; and that all degrees of saturation of the space considered 
are possible. 

The problem of dealing with such cases becomes tractable when we re- 
cognize that the behavior of each different fluid can be treated separately. 
Thus, for a specified fluid, there will exist at each point in space capable of 


being occupied by that fluid a macroscopic force intensity vector E, defined 
as the force per unit mass that would act upon a macroscopic element of the 
fluid if placed at that point. In addition, if the fluid does occupy the space, 
its macroscopic flow rate at the given point will be indicated by the velocity 
vector q, the volume of the fluid crossing unit area normal to the flow direction 
in unit time. 

We shall thus have for each fluid two superposed fields, a field of force 
and a field of flow, each independently determinable. The equations describing 
the properties of each of these fields, and their mutual interrelations, comprise 
the field equations of the system; these in turn, in conjunction with the boun- 
dary conditions, determine the nature of the flow. 


. The field of force-—We have seen already that the force per unit mass 
is given by 


E = g—(I/p)gradp, (15) 
and the force per unit volume by 

ta oF = eg — gradp. (17) 
; Either of these force vectors could be used, and either is determinable 
from the other, but before choosing one in preference to the other, let us first 
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consider the properties of their respective fields, of which the most important 
for present purposes is whether or not the field has a potential. To simplify 
our analysis we will make the approximations that 


Z = COnSt. (82) 


and for chemically homogeneous liquids under the range of temperatures 
and pressures normally encountered in the earth to drillable depths, 


© = const. (83) 
For gases, on the other hand, we shall have an equation of state 
e = f(p, T) (84) 


where T is the absolute temperature. 


The vector E has a particular value at each point in space and the ensemble 
of such values comprises its vector field. The criterion of whether this field 
has a potential, that is to say, of whether 


B= ora ©! 


where © is a scalar field, is whether the field E is irrotational, which can be 
determined from its curl. From equation (15), 


curl E = curl teas (1/¢) grad p] 
=V xe—yVv * [(1/p) vel. 
As is well known, the gravity field is irrotational even without the assump- 
tion that g = const, so that Vy x g =0. Also 
—V * [0/e) val = — yC/e) x ve — (1/9) YX We 


=— VU/p) x vp. 
Consequently 
curlE = — y(l/9) x yp = vp x VUl/e), (85) 
so that 
curl E = 0 when yp x y(l/p) = 0. (86) 


Therefore, in order for the field to be irrotational, and hence derivable from 
a scalar potential, it is necessary either that yp = 0, corresponding to constant 
pressure, or V(1/e) = 0, corresponding to constant density, or else that the 
vectors yp and (1/09) be colinear, corresponding to a coincidence of the sur- 
faces of equal density and equal pressure. 

The second of these three cases is satisfied by a liquid of constant density, 
and the third by a gas whose density is a function of the pressure only, such 
as occurs under either isothermal or adiabatic conditions. For the general 
case, however, of a gas for which 9 = f(p, T), and the surfaces of equal tem- 
perature do not coincide with those of equal pressure, then the surfaces 
© =const will also not coincide with the surfaces p = const and we shall 
have two intersecting families of surfaces, (1/e) = const, and p = const, 
for which equation (85) applies. 

This is the condition corresponding to thermal convection, and the fluid 
will have a convective circulation in the direction that will tend to bring the 
surfaces of equal density into coincidence with those of equal pressure, with 
the less dense fluid uppermost. 
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Fic. 8 — The potential ® as a line integral of the field of force E. 


Hence, subject to the condition that either 9 = const, or o = f(p), 
curl E-= 0 ‘and Ei == grad D- (87) 
The value of ® at any arbitrary point P in space (Fig. 8) is then obtained 
by é ) 
P 
@(P) = O(P,) — | Es ds 
Po 


A d 
= O(P,)—) | ge) | ds Ks (88) 
Po 0 a 


d Bd 
“ = QM(P,)+ez+ J sein 


Dp 
\ 
Po Po id 


= O(P,) + | gdz+ 


where the integral from P» to P is taken along any path s. Then by setting 
@(P,) = 0 when z = 0 and py = 1 atmosphere, we obtain 


P dp 
@(P) = ez 4+ | mp? (89) 


where p is now the gage pressure, or the absolute pressure less | atmosphere. 
If the fluid is incompressible and chemically homogeneous, this reduces 


to the simpler form 
@(P) = gz + pio. (90) 
For this case, if a manometer is tapped into the system at the point P, the 
height A above the level z = 0 to which the liquid will rise, will be 
h=z-+plog, (91) 
50 


from which it follows that 
gh =gz + pile =O, (92) 
in agreement with our earlier definition of ® in equation (14). 
If the fluid is incompressible and chemically inhomogeneous, as in the 
case of water of variable salinity, the density © will not be a function of 
pressure only, and, in general, surfaces of constant density will not be parallel 


to surfaces of constant pressure. For such a system curl E + 0, and no 
potential exists. 
We thus see that, with the exception of cases of thermal convection, 


and of inhomogeneous liquids of variable density, that the fields E for both 
liquids and gases are irrotational and are derivable from a potential ©. Since 


E is a force per unit of mass, then ® is an energy per unit of mass, and repre- 
sents the work required to transport the given fluid by a frictionless process 
along a prescribed (p,T)-path from a standard position and state to that 
_of the point considered. Surfaces ® = const are accordingly equipotential 
- surfaces, or surfaces of constant enérgy of position, and the fluid will tend 
to flow from higher to lower potentials or ‘energy levels. 


The field of force per unit volume, H, can be disposed of summarily. 
Since 


= eg — gradp, 
then 


curl H = Vv x (og) —V xX YP 
=VveXg+pyxe—V% w- 
But, since VY <x yp and Vv x g are each zero, then 
curl H = Ve x g. (93) 


This is.zero only when 9 is constant or when the surfaces of constant density 
are horizontal. The last condition never occurs except when the motion 
is vertical. Hence, for motion in any direction other than vertical, the field 


of the vector H does not have a potential except when the density of the 


~ fluid is constant. For the special case of constant density, 


ie gradIT , (94) 
where 
Il = p® = pgz +p (95) 
is the energy per unit volume of the fluid at any given point. 

In view of the fact that the field E has a potential for both liquids and 


' gases under the conditions specified above, whereas, in general, the field H 
‘has a potential only for the special case of liquids of constant density, then 
there is no advantage in using the latter in preference to the former, and 


henceforth it shall be dropped from further consideration. 


The generality of the field of force, as herein defined, merits attention. 


The force vector E for any given fiuid not only has values in space occupied 
. by that fluid, but also in any space capable of being occupied by the fluid. 


At a point in air, for example, the force Ew» for water would be 
= g—(1/Pw) grad p, 
of 


and since, in air, grad p is Pair g, then 
fe ft gee. (96) 
Ow 
The field E for a given fluid thus extends throughout all space of conti- 
nuous permeability. When several fluids are to be considered, then at each 
> 
point in space there will be a different value of E, for each separate fluid, given 
by: 


Ew = 


) 


E, —g—(l/o,)gradp, | 


En = g —(1/0n) grad p. | 


The vectors E for the separate fluids of different densities will differ 
among themselves, both in magnitude and direction, but will all fall in the 


same vertical plane, that defined by g and — grad p (Fig. 9). 
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FiG 9 — Force vectors E at the same point corresponding to fluids of different densities. 
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Similarly the potentials of different fluids at the same point will be: 
®, = ez +5 (@ie.), 


O,=gz4+) (dp/e2) , (98) 


®n = gz + J (dp/on), 
where (1, 2, ... Q, are the variable densities of the separate fluids. 
Since the equipotential surfaces for the separate fluids must be normal 


to the respective vectors E, then it follows that the equipotential surfaces 
of different fluids passing through a given point will not be parallel to one 
another, although they will all intersect along a line normal to the 


(eg. — grad p)-plane. 


The field of flow.—We have already defined the flow vector g for a space 
which is entirely filled with a single fluid. -For a space which is incompletely 
filled with a single fluid, or is occupied by two or more intermixed fluids, then 
there will be two or more superposed flow fields, not in general in the same 


— 
direction, and a separate value of q for each separate fluid. 

The principal condition which must Le satisfied by the field of flow inde- 
pendently of the field of force is that it must be in accord with the principle 
of the conservation of mass. Thus, if a closed surface S, fixed with respect 
to the porous solid, is inscribed within the field of flow, then the total net 
outward mass flux of any given fluid in unit time will te equal to the diminution 
of the mass of that fluid enclosed by S, assuming that processes which create 
or consume the given fluid are forbidden. 


Fic 10 — Mass discharge through a fixed closed-surface S as the basis for the equation 
of continuity. 
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This condition is expressed by 
§\ 9gn dS =—dm/dr, (99) 
Ss 


where g, is the outward-directed normal component of qg, and m the mass 
enclosed (Fig. 10). 

In the case of a space completely saturated with the given fluid, by 
dividing the integral (99) by the volume V and then letting V tend to zero, 
we obtain 

lim 1 y{ 1 om 
V= 0:04" aap nee 
which is the rate of loss of mass per unit macroscopic volume at a given point. 
Then, since 


div pq = pee. 


om/dt =fV Oop/odt), 
equation (100) becomes 
div og = —fQg/dr), (101) 


which is the so-called «equation of continuity» of the flow. If the motion 
is steady, then d¢/d¢ is zero, and the equation simplifies to 


div oq =0. (102) 
When this is expanded it becomes 
div eq = VV pg = Ve 9+ PY 9 =0," 
and when 0 is constant over space, corresponding to the flow of a homogeneous 
liquid, ye = 0, and 


divg =0. (103) 


Relation between field of force and field of flow.—In a space completely 
saturated by a single fluid, the field of flow and the field of force are linked 
together by Darcy’s law 

q.=oE, (104) 
which, in those cases for which curl B= 0, becomes 
q=—ograd®. (105) 


If the solid is isotropic with respect to permeability, the conductivity o 
is a scalar and the flowlines and the lines of force will coincide; if the solid 


is anisotropic, o will be a tensor and E and q will then differ somewhat in 
direction except when parallel to the principal axes of the tensor. 


Limiting our discussion to isotropic systems, by taking the curl of 
equation (105), we obtain 


curlg =—y xoVO =— Vo x VO —aV x VO. 


Then, since the last term to the right is zero, this becomes 


curlg = — yox yO, (106) 
which is zero only when c is constant throughout the field of flow. Therefore, 
in general, 

curlg £0, (107) 


and this circumstance precludes the derivation of the flow field from an assumed 
velocity potential, for, with the exception of the flow of a fluid of constant 


54 


density and viscosity in a space of constant permeability, no such function 
exists. 


The flow of a given fluid through a porous solid incompletely saturated 
with that fluid is equivalent to flow through a solid of reduced permeability, 
because the space available to. the flow diminishes as the saturation decreases. 
For saturations greater than some critical minimum value, the flow obeys 
Darcy’s law subject to the permeability having this reduced, or so-called, 
«relative-permeability» value. Thus, with two interspersed but immiscible 
fluids in the same macroscopic space, 


dQ SiGe, =~ or [g —(1/0,) grad p], | 
eer ates He ‘ (08) 
dz =Ore Ep = Gre[g — (1/p2) gard p], \ ‘ 


where o,, and o;, are the relative conductivities of the two fluids. 


It will be noted that, except for vertical motion, E, and E, are not parallel. 
Consequently the two force fields and flow fields will be, in general, transverse 
to one another. 

Pioneer work on relative permeability as a function of saturation was 
done on the single fluid, water, by L. A. Richards ® in 1931. Subsequently 
studies of the simultaneous flow of two or more fluids were initiated by 
Wyckoff and Botset, !° and by Hassler, Rice and Leeman * in 1936. Since 
that time many other such studies for the systems water-oil-gas have been 
published. 

One flaw which has been common to most of these multifluid experiments 
has been that the experimental arrangements and their interpretation were 
usually based upon the premise that the flowlines of the various components 
are all parallel and in the direction — grad p. Since this is far from true, 
there is some question of the degree of reliability of the results of such experi- 
ments. Even so, the existing evidence indicates that for saturations greater 
than. some critical minimum the flow obeys Darcy’s law in the form given 
in equations (108). For saturations less than this limit, there should still 


be a general drift of discontinuous fluid elements in the direction of the field E, 


but with, as yet, no well-defined relationship between q and E. 

The migration of petroleum and natural gas, through an otherwise water- 
saturated underground environment, from an initial state of high dispersion 
to final positions of concentration and entrapment, constitutes an example 
of the latter kind. 


Properties of the combined field of flow.—Al\though the present paper 
does not permit of their elaboration, practically everything that is now known 
concerning the flow of fluids in a porous-solid, three-dimensional space is 
deducible from the foregoing field equations. The equations for the flow 
of an incompressible fluid are of the same form as those for the steady conduc- 
tion of electricity. Consequently the well-known solutions to the electrical 
equations can be applied directly to analogous situations in fluid flow. The 
unsteady flow of a compressible fluid, when the force field is irrotational, 
is closely analogous to the unsteady conduction of heat, although somewhat 
more complex, and the solutions of the heat equations can be adapted with 
some modification to the analogous fluid-flow problems. 

The fluid phenomena which are unique and have no counterpart in other 
more familiar field theory are those involving multiple fluids. If we consider 
the flow of two immiscible fluids of unequal densities, such as water and_ oil, 
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interspersed in the same macroscopic space, we have seen from equation (97) 
—> 


that, in general, the force fields E, and E, of the two fluids will lie in the same 
vertical plane, but will have divergent directions, that for the less dense fluid 
being upward with respect to that for the more dense. As a consequence 
the two fluids will drift in these respective directions, and, in response to 
suitable impermeable barriers, will tend to become completely segregated. 

Once this segregation is achieved, for any further steady motion of either 
of the fluids, the interface will appear as an impermeable barrier. Across 
this interface, neglecting minor pressure differences due to capillarity, the 
pressure in the two fluids must be the same. Then, in case neither of the 
fluids is in motion, the interface will have to be horizontal, with the less dense 
fluid uppermost, since this is the only inclination of the surface for which 
the pressures on opposite sides can be the same. 

When either or both of the fluids is in motion, however, in a nonvertical 
direction, the vectors (1/9,) grad p and (1/p2) grad p will be inclined from 
the vertical, and the corresponding equipressure surfaces will be inclined from 
the horizontal by the same amounts. At the interface every equipressure 
surface in one system must match that of the same value in the other. 

At the same time, the flowlines in each system must be parallel to the 


interface so that the (E,, E,)-plane must be tangent to the interface. Conse- 
quently the interfacial surface and the flow patterns in the two systems must 
mutually adjust themselves until these conditions are simultaneously satisfied 
before a steady state of flow becomes possible. 

Of particular interest is the special case of this general situation wherein 
one of the two fluids is in motion and the other is completely static. In this 
case the equipressure surfaces in the static fluid are horizontal and surfaces 
equal difference of pressure are equally spaced; whereas, in the moving fluid 
the equipressure surfaces are inclined downward in the direction of the hori- 
zontal component of the flow. Their lines of intersection at the interface 
must accordingly be horizontal. Also, since the vector grad p lies in the 


same vertical plane as the vector i of the flowing fluid, it follows that the 
horizontal component of the flow direction must be parallel to the direction 
of steepest slope of the interface. 

If we consider a vertical plane parallel to the flow direction and perpen- 
dicular to the interface, then in this plane the surfaces of constant pressure 
will appear as lines of constant pressure, inclined in the flowing fluid and 
refracting into the horizontal across the interface in the static fluid. Let the 
flowing fluid be denoted by the subscript 1 and the static fluid by the subscript 2. 
Then along the interface in the direction of the flow 


(Op/ds), = (Op/ds)> (109) 
In the static fluid 
(Op/ds)2 = — Peg sin 9, (110) 


and consequently depends only upon the density 9, and the angle of slope 0, 
where @ is positive when the slope is upward in the flow direction. 
In the flowing fluid, since 


B= &s — (1/e)Op/ds), > 


- 


then 


(Op/ds)1 = Pigs — PiE 


=— igsin 0+ oe, |grad ®,|. { (111) 
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Equating (110) and (111) then gives 
ese sin 0 = 9,gsin § — p, |grad D,|, 


which, when solved for sin 9, becomes 


sin § = eee 
& P2— 1 
If the less dense fluid is flowing (Fig. 11a), 92 > 9,, the term to the right 
will be positive, and the interface will tilt upward; if the more dense fluid is 
flowing (Fig. 115), 91 > 0», the term to the right will be negative, and the 
interface will tilt downward. — 


\grad @,]. (112) 
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Fic. 11 — Interfaces between flowing and static fluids of unequal densities. 


(a) Upward tilt in the flow direction when less dense fluid is flowing. 
(6) Downward tilt when more dense fluid flows. 


Equation (112), derived earlier by Hubbert, 1% !8 is the fundamenta, 
equation pertaining to fresh -water — salt-water relations along shore lines, 
and to oil-water and oil-gas interfacial relations such as water or gas 
coning during oil and gas exploitation. It is also the basic equation 
governing the underground positions of oil and gas entrapment. If the water 
is static, accumulations of petroleum or of natural gas will occur beneath 


ayf/ 


downwardly concave impermeable barriers, with the oil-water or gas-water 
interface horizontal. If the water is in motion, as often is the case, the 
oil-water or gas-water interface will be inclined, and oil or gas may be trapped 
beneath structures completely unclosed in the hydrostatic sense. 


RESUME 


In paying our respects to M. Henry Darcy on this centennial occasion, 
we stated at the outset that it should be our endeavor: 

1. To show unequivocally what Darcy himself did and stated with respect 
to the relationship which now bears his name, and to give his results a more 
general, but still equivalent, physical formulation. 

2. To derive Darcy’s law directly from the fundamental Navier-Stokes 
equation of motion of viscous fluids. 

3. To develop, in at least their primitive forms, the principal field equations 
of the flow of fluids through porous solids. 

These objectives have now been accomplished, and the result is that, 
despite a number of troublesome complexities such as those arising from 
thermal convection and from waters of variable salinity, the field theory of 
the flow of underground fluids, is capable of being brought into the same 
kind of a comprehensive unification as that already achieved for the more 
familiar phenomena of electrical and thermal conduction. 

In closing it is pertinent to reiterate that Darcy’s empirical formulation : 


q=K(i—h)/l, 


which, as we have seen is valid for flow of a homogeneous liquid in any direc- 
tion, is physically equivalent to the expressions : 


q = —(Nd?\(p/p) g grad h, 


q=+ (Na?)(p/)1g — (1/¢) grad p] | 


From these it follows that fluids can and do flow in any direction whatever 
with respect to that of the pressure gradient. 

So far as we have seen no error of any kind with respect to the flow of 
fluids through porous solids can be attributed to Henry Darcy; and the errors 
which have been alleged by various authors during recent years appear on 
closer inspection to have been those committed by the authors themselves. 
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LES RECHERCHES DE TECHNIQUE PLUVIOMETRIQUE 
DE LA STATION CENTRALE SUISSE 
DE METEOROLOGIE, M. Z. A. ZURICH (complément) 


RESULTATS OBTENUS DU CHAMP D’ESSAIS 
DU SENTIER (JURA VAUDOIS), 1941-1951) 


par Paul-L. MERCANTON, Lausanne 
(Cf. P-V. de la session d’Oslo, aodt 1948: « Les recherches de technique 
pluviométrique », p. 131. Pour le détail voir: Annales de la M. Z. A., 1938 
1939, 1940, 1944 


Délimité au large des versants, sur le fond horizontal de la Vallée de 
Joux, en amont du lac et 4 1000 m d’altitude, ce terrain avait recu a fin avril 
1941, un ensemble d’engins pluviométriques étudiés céja a Zurich quant 
a leur neutralité aérodynamique, mais qu’il convenait de soumettre encore 
a des essais critiques en plein vent, dans des conditions variées et sévéres. 
A ce dernier égard, le climat de la Vallée, rigoureux en hiver, volontiers chaud 
a l’excés en été, assurait d’utiles expériences que la situation du champ, en 
travers des vents, humides, du secteur ouest a sud, garantissait pleinement 
aussi. 

Les observations se sont poursuivies sans discontinuer durant dix années 
par les soins de Mr. le professeur Pierre Baud principalement, avec la colla- 
boration du sussigné, pour la M. Z. A. L’équipement instrumental n’a guére 
changé durant cette décennie, laps de temps d’ailleurs nécessaire pour une 
comparaison pertinente des appareils. Les résultats ayant paru au complet 
dans les Annales de la M. Z. A., il suffira ici de les récapituler: 

Les deux grands totalisateurs du modeéle classique de la M. Z. A., (Maurer 
et Billwiller) occupant les extrémités de la diagonale du champ ont prouvé 
Puniforme ventilation du terrain: les quantités d’eau recueillies par eux et 
conservées d’automne a automne n’ont différé entre elles que de quelques 
dixiémes pour cent (moyenne des dix ans: 0,45 °%). Le modéle normal de 
la M.Z. A. a toujours recueilli un peu davantage d’eau que l’autre, dont 
l’écran était abaissé de 13cm par rapport a celui du premier dont l’écran 
a son bord supérieur au niveau de l’ouverture pluviométrique. Les essais au 
tunnel avaient avéré pour cet écran ainsi abaissé un optimum de neutralité 
dans le vent régulier, mais en est-il de méme dans le vent naturel? Dans le 
doute, et aucun des engins du champ ne pouvant étre tenu pour parfait, on a 
pris en général comme base de comparaison numérique la moyenne des indi- 
cations annuelles des deux totalisateurs. Bien entendu, on a réduit les indi- 
cations de tous les engins du champ a l’ouverture pluviométrique classique 
de 200 cm*. Notons enfin que le totalisateur a écran abaissé a emmagasiné 
en moyenne un peu moins de précipitations durant la saison neigeuse que le 
totalisateur normal, chez lequel I’écran, par sa position, favorise un peu len- 
trée de la neige. 

Le totalisateur ellipsoidal 4 vidange mensuelle, qui s*était montré aéro- 
dynamiquement parfait au tunnel, n’a pas répondu a notre espoir d’en faire 
la base des comparaisons avec les autres appareils. Si pour la pluie, le revé- 
tement de brosse métallique (I’« oursin »), a été certainement assez efficace, 
on n’a pu réaliser en hiver un chauffage rationnel de l’engin; le rejaillissement 
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des gouttes de la pluie ou le balayage des fiocons de la neige dans l’ouverture 
réceptrice ont en effet toujours exagéré les quantités d’eau recueillie. Nous 
avons donc da le disqualifier comme étalon. Il a mesuré en moyenne générale 
6 % de plus que les grands totalisateurs. Le pluvio-totalisateur mensuel M.Z.A.- 
Mercanton qui avait fait de bonnes preuves a Zurich, n’a pas supporté les 
rigueurs de l’hiver jurassien (—25° C et davantage) et a dd trés t6t étre supprimeé. 

Le Hellmann normal a mesuré 2 % en moyenne générale de plus que les 
totalisateurs annuels. Durant quelques années, mais en saison chaude seu- 
lement, on a expérimenté le pluviométre sphérique Haas-Liitschg, un globe 
métallique présentant au vent, quelle qu’en soit la direction, un ensemble de 
trous circulaires, symétriquement répartis, et faisant en projection sur le plan 
normal a la dite direction, la surface de réception de 200 cm”. L’eau engagée 
dans les ouvertures est canalisée vers un réservoir commun. L’appareil est 
affecté a un trés haut degré des défauts qui nous ont fait rebuter les obser- 
vations de l’engin ellipsoidal; il a toujours indiqué des hauteurs d’eau consi- 
dérablement supérieures a celles des pluviométres 4 ouverture horizontale, 
ce quise concoit d’emblée. La:signification des mesures au H-L n’est d’ailleurs 
pas claire. 

Un dernier dispositif mérite de retenir l’attention par ses résultats pro- 
metteurs: le pluviomeéetre Hellmann de haut modéle (60cm), entouré d’un 
écran de tdle plein centré sur l’axe du pluviométre et fait de deux troncs de 
cone circulaires, raccordés en fagon de poulie ou de diabolo par leurs petites 
sections. Haut de 60 cm aussi, cet écran présente en haut et en bas, des ouver- 
tures circulaires de 54 cm sans rebords et se rétrécissant graduellement jusqu’a 
44 cm en son milieu. Ces ouvertures terminales sont rigoureusement dans les 
plans de la bouche et du fond du pluviométre; l’ensemble est supporté par 
trois pieds équidistants. Les essais a la soufflerie ont démontré la neutralité 
excellente de cette construction. 

La comparaison pendant plusieurs années au champ d’essais avec |’Hell- 
mann nu a donné les résultats suivants: le pluviometre a écran en poulie a 
emmagasiné en pluie 2° de moins que |’autre; il a recueilli en neige 2% 
de plus. D’autre part ici encore une différence par excés est apparue avec la 
moyenne des totalisateurs annuels. Cette différence persistant entre les engins 
de mesures journaliéres du champ et ceux a longue échéance est troublante : 
elle force a se demander si la couche classique de 2mm d’huile, recouvrant 
la liqueur saline des totalisateurs, supprime bien toute évaporation? 

Les essais au champ donnent enfin un regain d’intérét a une question 
posée antérieurement par l’auteur (Recherches de technique pluviométrique 
III, Annales 1939): un pluviomeétre, abrité du vent par un entourage méme 
distant, d’arbres ou d’édifices élevés, ne recoit-il pas trop d’eau météorique? 
En effet, le pluviométre Hellmann normal du Collége installé a 200m de 
notre champ d’essais et 4 une quinzaine de métres au-dessus, a décelé pour les 
10 ans entiers, un excés de 3 % sur l’Hellmann du champ et de 5 % sur les 
totalisateurs annuels. 
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